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ABSTRACT 



A class of analytical and numerical solutions to the 
problem of inviscid-jet deflection from plane and axisym- 
metric concave surfaces are presented. The methods of solu- 
tion for the two-dimensional cases included a novel method 
developed in the present investigation, as well as those 
due to Levi-Civita and Riemann and Hilbert. These methods 
have provided indirect solutions for the shape of the curved 
solid boundary in terms of a given jet-departure and jet- 
deflection angle. The axisymmetric cases, namely, the deflec- 
tion of axisymmetric jets from hemispherical thrust reversers 
are solved directly through the use of the finite element 
method and a novel iteration scheme. The real-fluid effects 
on the jet deflection are estimated by comparing the results 
with those obtained experimentally. 

The analytical and numerical methods developed or used 
herein are sufficiently general to yield direct and indirect 
solutions for more general nozzle and curved-re verser combi- 
nations of special interest to thrust reversal on aircraft 
engines . 
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NOMENCLATURE 



A 

A e 

A i 

a 

a k 

a n 

B ? 

b 

b k 

C T 

a° t 

DELV 

d 

dA 

e 

FAC 

FACU 

FACL 

f 

f(t) 

H(t) 

h(a) 

I(4>) 

Im[ ] 



area of the meridional plane 

area of a triangular element e 

area of a subtriangle in a triangular element 

a parameter 

x j " x i 

coefficients of a series 
load matrix for element e 
a parameter 

r i " r j 

thrust coefficient 
error in C^ 

incremental correction to the free surface 
jet thickness 
elemental area 
denotes an element 
a multiplier 
a multiplier 
a multiplier 
a parameter 

shape function in Riemann-Hilbert method 
homogeneous solution of the Riemann problem 
function of a represented by a Fourier sine series 
functional of velocity potential 
functional I for an element e 
imaginary part of [ ] 
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i 


V-> 




elements of the stlfness matrix for element e 


k 


a parameter 


ill 
1 l* 2’ 3 


side lengths of a triangle 


M 


number of triangular elements 


M o 


a complex constant 


N 


number of terms In a series 


N o 


a complex constant 


P 


pressure 


P 

o 


pressure at the entrance of the nozzle 


Q(t) 


function of t in the Riemann-Hilbert method 


q 


magnitude of velocity 


R o 


radius of the lip of the bucket 


CD 

i — i 
i i 


real part of [ ] 


r 


radial coordinate, also modulus of £ 


r o 


radius of the nozzle 


r E » r F ’ r G 


radial coordinates of the points E, F and G 


r l> r 2> r 3 


radial coordinates of corner nodal points 1, 2 
and 3 in a triangular element 


S 


a complex variable 


s u 


an array related to the geometry of each element 


s 


separation between the nozzle and the bucket lip 


T 


thrust 




arrays expressing the geometrical characteristics 
of an element 


t 


a variable 


U,V 


horizontal and radial components of velocity 


e e 
u ,v 


velocities in an element e 
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V 

V 

o 

V c 

V J 

V 

s 

VC 

V0 

VS 

w 



x,y 



X q jY q 



a P s a q 
3 

C 



n 

n R 

0 

X 

P 

a 

4> e 

(11) a 
X i 



magnitude of velocity 

velocity at the entrance of the nozzle 
velocity at the point C 

final freestream velocity 3 (here V. = 1) 

velocity at a point over the assumed free streamline 

denotes V 

c 

denotes V 

o 

denotes V 

s 

complex potential function 
cartesian coordinate 
coordinates of the point Q 
a complex variable 

exterior angles of a polygon at P and Q 

deflection angle 

a complex variable 

area coordinates 

a dummy variable 

reverse-thrust efficiency 

direction of the velocity vector 

a constant 

density of fluid 

argument of x, 

velocity potential 

velocity potential in an element e 

specified normal derivative 

array of area coordinates 
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\p stream function 

a complex variable 

w a complex function of the velocity vector 
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I. INTRODUCTION 



The prediction of the deflection of finite jets by 
segmental or curved deflectors is a matter of importance 
in the design of targe t-type thrust reversers on jet-aircraft 
engines, of flip-buckets on spillways, of Pelton wheels, etc. 

In general, the deflection of a free jet by a solid boundary 
is well suited to potential-flow analysis because of the 
dominance of inertia and pressure intensity in the establish- 
ment of the flow pattern. The design of impulse machinery, 
thrust reversers, etc. utilizing this momentum change could 
be facilitated greatly if the idealized geometry of the system 
under potential flow conditions were known because such condi- 
tions represent asymptotic values which are approached as the 
effects of secondary variables such as entrainment, boundary 
layer, compressibility, jet attachment to adjacent surfaces, 
etc., are decreased. With such information available, refine- 
ments of design could be based upon a secure knowledge of the 
fundamentals, and many rules of thumb could be replaced with 
precise quantitative data. Specifically, if the total angle 
through which the jet is deflected is determined for conditions 
of both partial and complete interception by the boundary, then 
the principle of impulse and momentum can be used to compute 
forces or other dynamic characteristics of the system. 

The - problems discussed herein fall in the general category 
of finite or infinite, two-dimensional or three-dimensional. 
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viscous or inviscid, and compressible or incompressible flow 
past bodies of arbitrary shape. In fact, in its broadest 
sense, they constitute the most fundamental and also most 
complex problems of fluid dynamics. Cur inability to deter- 
mine the separation points, or lines, on a body and to estab- 
/ 

lish a link between the separation and the unsteady fluid 
motion in the wake of the body gave rise to numerous exact 
and approximate methods of analysis. Most of these analyses 
rely in part or whole on the experimentally observed charac- 
teristics of flow to determine the magnitude of the disposable 
parameters introduced into the analysis. Suffice it to say 
that future efforts will continue to be based upon partly a 
secure knowledge of the fundamentals of fluid mechanics and 
partly on experimental facts. 

Evidently, the introduction of simplifying assumptions such 
as inviscid, incompressible fluids, the use of bodies with 
a salient edge at which the separation point is fixed and 
finally the- restriction of the analysis to either two-dimen- 
sional or axisymmetric flows enable one to obtain solutions 
which are often of practical significance and in conformity 
with the experimental observations. 

The present work deals with inviscid, incompressible flows 
of finite extent impinging upon two-dimensional or axisymmetric 
bodies with sharp edges, and makes use of all the available 
exact or approximate methods of analysis in arriving at novel 
solutions towards the prediction of the characteristics of 
deflected jets. 
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A closer look at both the problem at hand and the methods 
available requires the further specification of the geometri- 
cal configurations of both the nozzle from which the jet is 
emerging and the body upon which the jet is impinging. The 
classical theory of jet deflection deals essentially with 
nozzles and bodies composed of straight segments. This is 
primarily because the successive transformations used in the 
analysis lead to relatively integrable equations for the 
segmental boundaries and to exceedingly complex integral 
equations in the case of curved boundaries. Furthermore, 
the fact that there cannot be a corresponding complex- function 
theory for three-dimensional space, the limitations of our 
understanding of the characteristics of axisymmetric or, in 
general, three-dimensional jet-deflection problems becomes 
quite clear. 

In order to bring into closer focus the results of the 
past efforts as well as the limitations of the existing 
methods, we will, in the foregoing, describe briefly the 
suggestions made concerning the "rounding of the corners" 
of solid boundaries, describe the essential features of the 
problems encountered through the use of a specific example, 
and then discuss the history of the two-dimensional jet 
deflection problem. The comments concerning the previous 
work on axisymmetric jet deflection will be reserved for 
Section III where such cases are studied through the use of 
the finite element method. 
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A. A BRIEF REVIEW OF SEGMENTAL JET DEFLECTORS IN 

TWO-DIMENSIONAL FLOWS 

The two-dimensional counterpart of the jet-deflection 
problem has been treated by several investigators through the 
use of the powerful analytic-function theory and successive 
conformal transformations. Sarpkaya [1] solved the U-shaped, 
two-segment, deflector problem where the turning angle between 
the segments is limited to SO degrees. Tinney et al. [2] 
extended this analysis to the case where the turning angle 
between the symmetrically situated segments is greater than 
90 degrees. Later, Chang and Conley [ 3 ] presented an analysis 
for a bucket composed of a series of segments of arbitrary 
number, length, and angles; however, the basic as well as 
practical problem of the direct analysis of jet deflection 
by curved buckets remains unsolved primarily because of the 
mathematical difficulties encountered. It is partly for the 
purpose of illustrating these difficulties and partly for 
the purpose of familiarizing the reader with the fundamental 
concepts and the analytical transformations that the basic 
steps in the analysis of the deflection of a two-dimensional 
jet from a U-shaped segmental deflector will be discussed. 

The flow in the physical z-plane (see Fig. 1) may be 
transformed into a hodograph-plane through the use of 



e 



-u + iv 




dW 

dz 



( 1 ) 



The flow in the hodograph plane may be transformed by means 






(a) z-plane 



(b) f -plane 



(c) fi-plane 




(d) t-plane 



•1 -k 0 k 1 ft 

r\ r\ r\ s?\ : 

D E F ABC 



Figure 1 Transformation planes for the U-shaped bucket 
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of Planck's logarithmic transformation into an ft-plane 
where ft is defined by, 

V. 

ft = ln(- -J-) = - In + i0 (2) 

* J 

The ft-plane is composed of straight lines. The horizontal 
lines represent the solid boundaries and the vertical lines 
the free surfaces. Such a polygon may be transformed (e.g. 
see Milne-Tomson [4]) through the use of the Schwarz-Chris- 
toffel transformation onto either the upper or the lower 
half of a t-plane by writing 



ft = M 0 / 

(t-t 1 ) 



dt 

oTTtt cxq/tt 

1 (t-t 2 ) ^ . 



+ N„ 



(3) 



Three of the t values in the t-plane may be chosen arbitrarily 
in accordance with the Schwarz-Chris tof fel transformation and 
the rest must be assigned unknown parameters. Accordingly, 
the points E, A, C are assigned the values -k, k and f 
respectively. Then, Eq. (3) becomes 



ft = Mo / — 5 + No 

(t -k'j V t -1 



(4) 



After integration and evaluation of the constants M 0 and 1^, 
one obtains 



ft = In 



Vw? 



+ V i-t 2 



y] t^-i 



for t _> 1 and t £ -k 

(5) 
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and 



ft = In 



Vl-k 2 + V 1-t 2 

-.\TW 



for -k < t < 1 



( 6 ) 



The potential function W in the upper half of the t-plane 
may be written in terms of the sources and sinks representing 
the flow. Noting that there is a source at D and a sink at 
C and that only the upper half of the t-plane is used, one 
has 

V.d 

W = -*2— [in (t-f) - In (t+1)] (7) 

IT 



From this relationship and the definition of ? in Eq. (1) 
and ft in Eq. (2) one obtains the final expression that relates 
the z-plane to the t-plane, 

z _ f \I 1-t 2 + V 1-k 2 r 1 It.. , 0 , 

a ' / r~o- f [ FT • t+I ] dt (8> 

J V k 2 -t 2 

N 

Since the purpose of this section was only to show the general 
procedure used for this type of transformations, the results 
of the integration are not shown here. The details are given 
by Sarpkaya [1]. 

It is apparent from the foregoing that whenever the solid 
boundaries in the physical plane (see Fig. 1) are composed of 
straight segments, then the ft-plane also is composed of 
straight lines. In other words, the flow retains a constant 
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direction along the solid boundaries (lines parallel to the 
logarithmic axis in the fi-plane), and a constant velocity 
along the free surfaces (lines normal to the logarithmic 
axis in the P-plane). Since such a polygonal boundary can 
always be transformed onto a t-plane through the use of the 
Schwarz-Christof fel transformation, all two-dimensional jet- 
deflection problems of this nature may be, at least theoreti- 
cally, solved. 

It is also apparent from the foregoing that whenever both 
the magnitude and the direction of the velocity vary along a 
rigid boundary, then the ft-plane is not in general composed 
of straight lines or of lines conveniently intersecting each 
other at 90 or 180 degree angle. In general if the curved 
boundaries in the physical plane is given, then the corres- 
ponding portion of the Q-plane will be an unknown curve 
since the magnitude of the velocity is not initially known. 
One may, therefore, attempt to simplify the problem, for 
example, by specifying the variation of the magnitude of the 
velocity along the boundary. This does not in any way alle- 
viate the difficulties to be encountered later since there 
does not exist a transformation function comparable to that 
given by Schwarz and Christoffel which could transform an 
ft-plane composed in part by curved segments onto a t-plane 
in which the potential-function may be written. Consequently, 
one will either seek other methods of handling the curved 
boundaries or completely abandon the direct approach of 
obtaining a solution for a given geometry. The indirect 
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approach called into action will seek a family of solutions 
for a family of initially unknown curved boundaries and then 
will let the designer choose one among those obtained to suit 
his desired needs. In the following, the efforts made by 
others in rounding the corners of buckets otherwise composed 
of straight lines and the efforts to obtain indirect solutions 
will be described. 

B. ROUNDING OF THE CORNERS AND BOUNDARIES 

Because of its theoretical as well as technological im- 
portance, the problem of rounding the sharp corners on rigid 
boundaries over or through which the fluid flows has attracted 
the attention of many mathematicians and hydrodynamicists . 
Prior to a brief exposition of these efforts, it should be 
noted that there has not been, up to now, a sufficiently 
general method of studying the flow around curved bodies and 
that the techniques so far developed are not only difficult 
to handle but also require the use of special techniques for 
a given problem. ■ 

1. Curved Walls of Constant Pressure 

As cited earlier, both the magnitude and the direc- 
tion of velocity vary along a curved boundary and the func- 
tional relationship between the wall geometry and the vector 
velocity is not known a priori. However, by making a simpli- 
fying assumption such that the pressure or the total velocity 
along thd curved wall remains constant, one can round the 

i 

corners in a bucket otherwise composed of straight segments. 
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The method may best be described by applying it to the U- 
shaped bucket case discussed earlier. 



Let the ft-plane (see Fig. lc) be modified in such a 



manner so that it takes the shape shown in Fig. 2a. Evidently 
the magnitude of velocity along PQ remains constant and the 
direction varies from tt/2 at P to tt at Q. The ft-plane is 
still composed of straight lines and may, therefore, be 
transformed into either the upper or the lower half of a 
t-plane through the use of the Schwarz-Christof fel transfor- 
mation. The remainder of the analysis will not be presented 
here. Suffice it to say that the transformation now involves 



therefore, renders the succeeding integrations relatively 
more difficult to perform. Be that as it may, the method 
is in principle capable of generating curved corners along 
which the pressure or the velocity remains constant. There 
is sufficient flexibility in the analysis to move the points 
A and Q along the boundaries to obtain smaller or larger 
curved corners within reasonable limits. 

2 . Curved Walls with Special Pressure Distributions 



is its modification in such a manner that while the direction 
of flow along the curved wall varies linearly, the magnitude 
of the velocity varies logarithmically. This procedure still 
allows the ^-plane to be composed of straight segments (see 
Fig. 2b). It is immediately apparent that the application 



two additional terms in the form of 




and 



A bolder approach to the modification of the ft-plane 
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B 

C 




F 




(a) ft -plane for 

constant pressure 
curvature 



E 



— ■ 1 

-ln(q/Vj) 



(b) 



B 

C 



D 

ft-plane for 
variable pressure 
curvature 





(c) ft- or z-plane (d) t-plane with modifications at P and C 



Figure 2 Modifications of the ft- and t-planes to round the corne: 
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of the Schwarz-Christof fel transformation now yields terms 
such as 



a /it a /it 

(t-t ) p (t-t ) q 



as a consequence of the fact that the line segments in the 
ft-plane, no longer intersect each other normally. Only 
under rare circumstances, i.e. for special values of 
or otq, that the integral equation resulting from the Schwarz- 
Christoffel transformation may be integrated in terms of 
the Beta-functions. In general, the transformation may be 
integrated only numerically. Since the determination of 
the coordinates of the rigid boundaries in the physical 
plane requires another integration, similar to Eq . (8), the 
solution of the problem via a two-stage numerical integration 
becomes not only difficult but also quite approximate. 

This method has been applied by Schieldrop [5] to 
the analysis of Borda-mouth-piece type flows with curved 
edges. The method is not in general easy to apply and requires 
a great deal of ingenuity and sophistication in numerical 
integrations . 

3. Curved Walls Through the Modification of the t-plane 

Cockcroft [6] introduced a method by means of which 
the corners of a solid boundary may be slightly rounded by 
modifying the t-plane. Consider the flow in a 90 degree 
corner as shown in Fig. 2c. This plane may be regarded as 
either z or the ft-plane. In the absence of round corners, 
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the Schwarz-Christof fel transformation relates this plane to 
the t-plane (see Fig. 2d) by 



z or p 



= Mo /* J- 
J (t- 



t — 0 dt 



(t-a) Y t-b 



+ N n 



(9) 



The corner at B may be rounded in an indirect way by replacing 
J t-0 in the above equation by \j t+1 + X \J t-1 . In this 

substitution the argument of \j t which changed to tt/2 at 
t = 0 now varies smoothly when t moves from t = -1 to t = +1. 
In other words, a curve is produced in the z or ft-plane whose 
curvature depends on the arbitrary value of X. It should be 
noted that for a value of X, say 



X 




the velocity or the pressure along the curved wall remains 
constant. In fact, evaluating dW/dz in the range of -l<t<l, 
one obtains 




This method was successfully used by Cockcroft [6] 
in calculating the electrical-stress distribution around 
certain two-dimensional conductors. A brief discussion of 
this method and the difficulties associated with its appli- 
cations are also described by Carrier et al. C 7 3 • 
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*1 . Curving of the Boundaries Through the Use of 
Special Transformations 

The Riemann-Hilbert transformation has been exten- 
sively used by Larock and Street [8] in the analysis of 
flow about supercavitating hydrofoils and by Larock [9] in 
analyzing in an indirect manner the efflux from curved 
nozzles. The mathematical details of this method will be 
discussed further later in connection with its application 
to the jet deflection problem. 

5 . Analysis of Flow about Curved Bodies by 
Levi-Civita 1 s Method 

Because of its challenging nature, the question of 
the existence and uniqueness of potential flows past general 
obstacles, having boundaries of given shape, has intrigued 
many outstanding mathematicians such as Brillouin [10], 

Villat [11], Levi-Civita [12], and others. In fact, Levi- 
Civita may be said to have solved the inverse problem of 
describing the class of all jets divided by curved barriers. 
Although no specific solutions have been presented, the prob- 
lem of determining such flows was reduced in special cases 
to the solution of nonlinear integral equations with appro- 
priate boundary conditions. The solutions are obtained in 
general by solving such equations numerically and by suitably 
choosing the necessary number of parameters. Such parameters 
described for example the variation of the inclination of 
the velocity vector, the variation of the magnitude of veloc- 
ity, etc. It is important to know that the present problem 
of the deflection of a two-dimensional jet from a curved 
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boundary has not been previously solved by any of the methods 
so far described. 

In Section II, Levi-Civita* s method, among others, 
will be described in greater detail in generating a family 
of curved buckets for a given jet and prescribed jet 
deflection angle. 

C. SUMMARY 

Major emphasis in the foregoing introduction has been 
on the direct or indirect analysis of the impingement of 
ideal fluids on two-dimensional straight or curved boundaries. 
These analyses are often difficult to apply and the amount 
of work required may not be commensurate with the need for 
a solution. In such cases, the use of relatively less- 
exact methods may be adequate or even necessary. For two- 
dimensional flows in general and for axisymmetric flows in 
particular approximate methods have been developed. A de- 
tailed discussion of these methods with the exception of the 
finite element method may be found in Robertson [13]. These 
methods may be classified as graphical (field plot or flow 
net), finite-difference numerical analysis (iteration, relax- 
ation, Monte Carlo), and analogies (electrical analogy, Hele 
Shaw, and membrane analogy), and finally the finite element 
method. The first three of these approximate methods have 
been widely used and amply described by Robertson [13]. The 
finite element method which came into existence during the 
past decade is currently being applied to the solution of 
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two-dimensional and axisymmetric Laplace-field problems. 

A detailed discussion of this method and its application to 
axisymmetric jet deflection problems will be presented in 
Section III. In the following, exact methods of analysis 
of the deflection of two-dimensional jets from curved buckets 
will be undertaken. 
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II. ANALYSIS OP TWO-DIMENSIONAL JET 



DEFLECTION FROM CURVED BUCKETS 

In the following, three methods will be employed to 
obtain a family of curved-bucket shapes in terms of a given 
jet deflection angle and the angle of departure at the lip of 
the two-dimensional bucket. The jet will be assumed to be 
free and impinging symmetrically upon the bucket. These 
methods consist of the modified hodograph-method, Levi- 
Civita's method, and the Riemann-Hilbert method. The first 
of these is novel and the other two have been previously 
employed as discussed in Section I.B. 

A. MODIFIED HODOGRAPH METHOD 

The method consists of the definition of the hodograph- 
plane, its transformation into an infinite strip, and subse- 
quently, into a t-plane through the use of the Schwarz-Chris- 
toffel transformation. The flow in the t-plane is expressed 
in terms of appropriate singularities and then the geometry 
of the curved bucket is determined through the evaluation 
of the proper integrals.. 

Consider the hodograph-plane in Fig. 3b where 



_ 1 dW q -i6 

5 = - vT a? = v: e 

J J 



( 10 ) 



The fluid stagnates at A, accelerates along ANB in a manner 
specified by the hodograph-plane, leaves the bucket at B at 
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iS 2 


B 


c 


M 


D 




M' 


C' 


B’ 



















(c) S -plane 



- Vk 

B ' C ' 



D 



- k 



0 

■^fc 



C B 



dx 

A 



t 






(d) t-plane 



Figure 3 Transformation planes for the modif ied-hodograph method 
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an angle of it, and eventually acquires an angle of deflection 
3 . Thus, the outer unit circle in the hodograph-plane repre- 
sents the free streamlines of constant velocity. As stated 
earlier, the fluid is assumed to be inviscid and the real- 
fluid effects such as boundary layer, entrainment, etc. are 
ignored . 

Clearly, the arc ANB in the hodograph plane may be chosen 
to represent different velocity variations along the bucket 
and different angles of departure at B, e.g., the curve AE. 
Although the analysis will be relatively more complex, there 
are, to be sure, no conceptual difficulties and differences 
between that presented herein and other cases. 

The hodograph-plane may be transformed to an S-plane 
(see Fig. 3c) through 



S 



= _ 2 i C 
C + 1 



(ID 



and the interior of the polygon in the S-plane may be trans- 
formed into the lower half of the t-plane by 



S = — In t 

7T 



( 12 ) 



which is obtained through the use of the Schwarz-Christof fel 
transformation. 

The flow in the t-plane is comprised of two sinks and a 
source. Hence the complex potential is given by 
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( 13 ) 



V.d , 

W = — [ln(t+k) + ln(t+^) - 2 ln(t+l)] 



Writing 



dz = - i - dW 

V j 5 



to obtain the geometry of the bucket, one has 



(jz = — r 1 + I • r_ — + —2: _ ? 

a tt L1 In t J L t+k , .1 t+1 

t+ k 



] dt 



( 1*0 



After separation of the real and imaginary part, one obtains 



x 1 , (t+k) 

d = ? ln — a 



(t+ir 



x 



( 15 ) 



and 



2. , 
d 



-2 / 
1 



ln t t+k 



<4^ + 1 



t+ ir 



t+1 



) dt 



( 16 ) 



The deflection angle 3 is related to the parameter k through 



k = EXP (- it -- 1 - n - £ ) 
1 - cos 3 



( 17 ) 



““ i R 

which is obtained by writing z, = e , inserting it in 
Eq. (11) and combining with Eq. (12). 
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Figure 4 Family of buckets obtained via modif ied-hodograph method 



The foregoing completes the formal analysis of the problem. 
The bucket geometry may now be calculated for a given deflec- 
tion angle & through the use of Eqs. (15), ( 16 ) and (17). 

The results obtained through a s traightf orward integration 
scheme are shown in Fig. k for a set of representative jet 
deflection angles. The coordinates of the 50 degree case 
are shown in Appendix A, and the computer program in Appendix 
B. Suffice it to say that the thrust acting on the bucket 
may easily be determined through the use of the momentum 
equation in terms of the deflection angle 0. 

B. LEVI-CIVITA 1 S METHOD 

This method, which is fully described in several reference 
texts (e.g. Birkhoff and Zarantonello [14], Milne-Thomson 
[4]), consists of the definition of a complex function w as 

w = 0 + i In (^-) 

J 

where 6 and q are’ respectively the direction and the magnitude 

of the velocity vector. V. is the constant velocity along 

J 

the free streamlines. The w-function is expressed in terms 
of a polynomial in such a manner that w is real over the free 
streamlines and complex over the curved obstacle. The shape 
or the curvature of the obstacle is defined in terms of the 
coefficients of a polynomial representing w. 

Even though this method was originally devised to deal 
with curved boundaries, as discussed in Section I.B-5, its 
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application has been mostly restricted to boundaries with 
straight-line segments (Cisotti [15]) , because of the ana- 
lytical difficulties encountered in its application to curved 
boundaries. Wu [l6] has employed it in the study of flow 
about fully-cavitating, slightly-curved, hydrofoils. 

The physical z-plane (see Fig. 5a) consists of a jet of 

width 2d and velocity V. impinging upon a curved, symmetric, 

0 

two-dimensional bucket BAB' . The complex-potential function 
W is given as usual by (see Fig. 5b) 

W = 4> + i^ (15) 

A s traight forward application of the Schwarz-Christof fel 
transformation establishes a relation between the t- and 
W-planes . Thus, one has 

V.d 

W = - — “ — [In (1-t) + In (1+t)] (19) 

TT 

Because of the symmetry, t is taken equal to +1 and -1 at 
C and C' respectively and zero at A. The parameter k is 
assigned to B and, as will be shown later, it is uniquely 
related to the deflection angle 3. 

The upper half of the t-plane is transformed into the 
inner region of a semicircle of unit radius (Fig. 5c) through 
the transformation given by Levi-Civita [12] as 

t = - | (C + |) (20) 
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(c) t-plane 



(d) f-plane 



Figure 5 Transformation planes for Levi -Civita ' s method 



3 ^ 



The circular arc BAB' in Pig. 5c represents the solid boun- 
dary over which the direction of the velocity vector is 
assumed to be prescribed and diametric line BCDC'B' repre- 
sents the free streamlines where the magnitude of the velocity 
is known. 

The relationship between the points in the z, t and 
C-planes are given by: 



point B 1 , t = -k, 

point A , t = 0, 

point D , t 

point B j t = k. 



C = +1 
C = i 
C = 0 
C = -i 



( 21 ) 



The W-plane is related to the z-plane, as before, by 




dW _ -i6 

dz V. 6 
J 



( 22 ) 



where q and 6 represent respectively the magnitude and direc- 
tion of the velocity vector. 

One can now define an w-function which is quite similar 
to that given by Planck as 



w = i In (-^- ^) = 0 + i In (^-) (23) 

O' j 

Nov-/ solving for z in the above equation and expressing all 
the variables as functions of c one has 



35 



z_ 

d 




iw( c) 




( 2 * 1 ) 



Defining 






r e 



ia 



and noting that over the bucket only the argument and along 
the free streamlines only the modulus varies, the combination 
of Eqs . (19) and (24) yields. 



z_ 

d 



1 

IT 



f e iw(a) 




cosa sina 

,2 2 
k cos a 



] da 



(25) 



and 



z 

d 



1 

IT 



2 2: 

/e lM(r) [ ? i- (r ~ r 2) . 



-] dr 



(26) 



^-(r + i ) 2 - 



~*r 



Equation (25) gives the coordinates of the bucket and Eq. (26) 
the coordinates of the free streamlines. 

The function co is assumed, according to Levi-Civita, to 
be of the form 



w(c) = i In ^ + ^ + a Q + a ][ c + a 2 ? 2 + a^C^ •••• ( 2 7) 

where the logarithmic term can be identified as that corres- 
ponding to a free jet impinging normally upon a flat plate 
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(Cisotti[17] ) • The additional terms in the polynomial whose 
coefficients are to be prescribed later modifies the expres- 
sion for a) and helps to curve the solid boundary. 

Because of the symmetry of the flow in the physical 
plane, the coefficients a Q , a zj> •••• are not acceptable. 

Noting that over the solid boundary r = 1 and over the 
free streamlines a = 0 or a = it, Eq. (27) may be reduced 
over the bucket to 



w 



(a) = i In 



| cos a | 

1 + sine 



*f + 



ia . i3o . i5c 
a,e + a 0 e J + a,-e 
lo o 



( 28 ) 



where plus sign is to be used for ir > a > tt/2, and minus 
sign for tt/2 > a > 0. Over the free streamlines, it reduces 
to 



w(r) = - tan 1 — — -5 + a_r + a 0 r^ + a c r^ ... ( 29 ) 

1 - r 2 1 3 5 

The coefficients a^, a^, a^, ... may now be chosen to obtain 
different bucket shapes. In the present study, w is chosen 
to satisfy the following conditions 



w( + l) = -it 
w(i) = i | - i ® 
w(-l) = TT 
to(o) = 0 



(30) 
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These conditions simply state that the jet departure angle 
at the lip of the bucket be equal to ±tt, and that the bucket 
be continuous in curvature at the stagnation point A and 
symmetric with respect to the jet axis. 

There are several ways to select the coefficients a^ . 

In this study a Pourier-series representation was used 
primarily because of the versatility afforded by it in 
representing different bucket shapes. 

Equation (28) can be written over the circular arc AB 1 
as 



w(a) = i{ln t — gos . - g + a, sino + a_ sin 3a + a. sin 5a...} 
J- t sma JL j 

(3D 

+{- j + a 1 cos a + a^ cos 3a + cos 5a ...} 

The real part of Eq. (31) represents the direction 0 of the 
velocity vector over the bucket. Thus, it is obvious that, 
one can prescribe the curvature of the bucket or the modulus 
of the velocity vector by properly choosing the coefficients 
a^ through a . 

Let 0 be represented by the lines shown in Fig. 6. Evi- 
dently, by choosing sufficient number of terms, one can re- 
present the variation of 0 with a shown in Fig. 6. In doing 
so it is important to make 0 symmetric with respect to the 
0 axis since the series are expressed in cosine terms only. 

Rewriting Eq. (31) one has over the arc AB* 
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Figure 6 Approximate representation of 0(a) 
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w(a) = i In + i h (a) + 0(a) 

1 t sin a 



( 32 ) 



where 



N 

h(a ) = Z a sin no (33) 

n=l n 



and 



0(a) 



N 

tt + E a cos na 
2 , n 

n=l 



(34) 



whe re 



a = 
n 



(o x - a 2 )n 



2 [cos na^ - cos na 2 ] 



(35) 



in which a 1 and a 2 represent the two parameters shown in 

Fig. 6. The evaluation of the Fourier coefficients a is 

shown in Appendix C. To complete the analysis one needs to 

iw(a) 



evaluate the term e 
of Eq. (32). Writing 



which appears in Eq. (25) in terms 



e iw(a) _ 1 + sin a e ~h(a) i0(a) 
cos a 



and separating the real and imaginary parts, one has 



ico(a) 1 + sm a -h(a) , N 

; = e (cos 0(a) + i sin 0(a)) (36) 

cos a 
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Now combining Eqs . ( 25 ), (33), (35) and (36), one finally 
obtains 

| = i f° e- h(o) sin e(o) 2k2 s . ln % sln 0) do 

d ’ j/2 1 - k 2 cos 2 o 

(37) 



and 



2.-1 
d it 



tt/2 



/ e- h(0) cos e(o) . 2k - l .;n g ( , l ._ t si . n _ o) 



2 2 
1 - k cos o 



do 

(38) 



Equations (37) and (38) yield the coordinates of each 
segment of the curved boundary through the use of the appro- 
priate values of o-^ and . This in essence completes the 
application of the Levi-Civita ’ s method to the analysis of 
the jet deflection from curved obstacles. In the following 
the implications of the variations of o^ and will be 
discussed . 

The variation of 0 with a shown in Fig. 6 may be written 
for each interval as 



0(o) = ~tt 0 < 0 < 0^ 

0(a) = \ ( 2 ) - f (2 + — ^ ) o 1 < 0 < a 2 

0 2 _ 0 ^ °2 ~ a l 

0(o) ~ ~ \ 02 <a< Tj (39) 

0(o) ~ ~ \ ^-<0<rr-0 2 



Hi 



0 ( 0 ) 



tt / a 
2 - a. 




T\-a^<a<T\-a^ 



0(a) = 0 



it - < a 



< n 



It will be shown that this particular functional depen- 
dence produces bucket geometries of the form shown in Fig. 7» 
Rewriting Eqs . (37) and (38) for the particular example 
under consideration one has 



x 

■f-o 



(40a) 



and 



h = % s 2 e - h(o) r2 sln g a * |ia CT) d 0 (not,) 

a 17 tt/ 2 1 - k* cos^a 



for the interval of ^ > a >o^> i.e. for the segment AP } 
(see Fig. 7) 



2. = 2. / e - h(o) k2 sln °< 1 4 gln g) oos[| da 

a 2 1 - k cos a 21 



(4la) 



and 



y 

d 



1 ,V h(o) k 2 slnatl + sir, a) sln rt ° ~ °1 .] dr + y 
* o 2 1 - k 2 cos 2 a 2 °2 - °1 p 

(41b) 
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Figure 7 Schematic representation of three different bucket shapes 
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Figure 8 Family of curved buckets obtained by Levi -Civita ' s method 



i.e. for the curved segment 



for the interval o 



2 



> a > o 



1 



PQ, and finally 




-h(a) sin (1 + sin a) 
2 2 

1 - k cos a 



do + x 



q 



(42a) 



and 




(42b) 



for the interval o^ > o > 0, i.e. for the segment QB. 

The coordinates of the free streamlines may be obtained 
by integrating Eq . (26) through the use of Eq. (29) and the 
coefficients a^ given by Eq. (37)- Finally, the deflection 
angle is evaluated by replacing in Eq . (29) the corresponding 



The foregoing equations have been numerically integrated 
(the computer program is given in Appendix D) , and the re- 
sulting family of 'bucket shapes are shown in Fig. 8a. A 
special case for a = 0.31416, o^ = 0.47124, and k = 0.946 
is shown in Fig. 8b. 

C. RIEMANN-HILBERT METHOD 

The method consists of the transformation of the physical 
plane into the W-plane, where W is given by 



value of r at the point C 




). 



W = <j> + i 
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and the transformation of the W-plane into the upper half 
of a t-plane through the use of the Schwarz-Chris tof fel 
transformation. Then the flow in the t-plane, which is 
composed of segments of solid boundaries and free stream- 
lines are related to the Riemann-Hilbert transformation as a 
well-posed mixed boundary value problem. Then a certain 
shape function f(t) is 'assumed and the characteristics of 
the flow in the physical plane are evaluated through 
integration. 

The details of the technique will become clear as it 
is applied to the determination of the deflection angle of 
a two-dimensional jet impinging symmetrically upon a 
two-dimensional bucket. 

Figure 9a shows a jet of width 2d impinging upon a curved 
bucket with a 180 degree angle of departure at the upper and 
lower lip of the bucket. 

It is apparent from the application of Bernoulli's equa- 
tion that the velocity is constant along the free surfaces. 

The complex- function W, for the upper half of the physical 
plane, is represented by an infinite slit as shown in Fig. 9b. 
The W-plane can in turn be transformed into the upper half 
of a t-plane through the use of the Schwarz-Chris tof fel 
transformation. The s traightforward application of this 
transformation yields 



W(t) 






IT 



In t 
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(42) 



Figure 9 
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Transformation planes for the Riemann-Hilbert method 
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where is the velocity along the free surfaces. In 
accordance with the rule of this transformation, three of 
the t values are assigned arbitrary values (t = 1 at A, 
t = +“ at C and t = 0 at D), and an unknown parameter k is 
assigned to the point B. 

The complex-potential-function W is related to the 
physical z-plane by 



dW 

dz 



-q 




(43) 



where q and 0 represent respectively the magnitude and the 
direction of the velocity vector. 

It is convenient at this point to separate the magnitude 
and the direction of the velocity through the application 
of Planck's logarithmic transformation given by 



w = In - i 0 (44) 

J 

Combining Eqs. (42) and (44) and solving for z in terms of 
t, one has 



f = - i / e- w(t) (45) 

It is apparent from the foregoing development that the real 

part of the complex-function w(t) is zero along the free 

streamlines where q = V.. Thus, 

3 



Re[w(t)] = 0 for t > k and t < 0 



The imaginary part of co(t) which represents the direction 
of the velocity vector is zero along the line AD where 
0 < t < 1. Over the bucket, however, both the magnitude 
and the direction of the velocity are variable and are an 
unknown function of t‘‘. Leu this function be denoted by f(t) 
in the interval 1 < t < k. 

The foregoing conditions define a well posed Riemann- 
Hilbert mixed-boundary value problem where the real and 
imaginary parts of an analytic function w(t) are alternately 
known along the real axis, i.e., 

Re(w)=0 t > k 

, Im(to) = -f(t) 1 < t < k (46) 

Im(w) = 0 0 < t < 1 

Re(io)=0 t < 0 

The general solution of the Riemann-Hilbert problem is well 

known (Cheng and Rott [18] and Larock [9]). 

In principle one writes 

w(t)=H(t)Q(t) (47) 
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where H(t) is the solution of the homogeneous-Riemann prob- 
lem (i.e. for f(t) = 0) as it satisfies all the other condi- 
tions listed in Eq. (46) and may be chosen as 



H(t) 




(48) 



In passing it should be noted that the particular form of 
H(t) which constitutes a solution for f(t) = 0 is not of 
special importance as shown by Song [19] and that uniqueness 
of the solution does not depend on the particular value of 
H(t ) . 

The function Q(t) may be written as (Larock [9]) 

Q(t) = - / dn (49) 

TT X) - t 

nn 1 



Combining Eqs. 



w(t) 



(46), (47) and (48), one has 

V -t( t - k) * -f (n) 

11 1 -n(n - k) 



dn 

(n - t) 

(50) 



The remainder of the analysis depends on the selection 
of f(t) . As pointed out earlier, f(t) would have been known 
had the shape of the solid boundary been known a priori. The 
present analysis, however, is an indirect approach and seeks 
the shape of the boundary in terms of the boundary conditions 
assigned and the deflection angle assumed. Thus, it is 
apparent that one may select f(t) as an arbitrary function 
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of t provided that the boundary conditions at the stagnation 
point and at the lip of the bucket are satisfied, and obtain 
a family of bucket shapes. In the following, f(t) will be 
taken as a linear function of t partly for the sake of sim- 
plicity and partly in order not to obscure the application 
of the method by introducing rather complex and unmanageable 
integrals. Thus, writing, 



and choosing a and b in such a manner that the slope of the 
bucket be ± tt/2 at the point A and ± it at the lip, one has 



f(t) = a + bt 



(5D 



it k - 2 
a 2 k - 1 



and 




Inserting Eq. (51) in Eq. (50) one obtains 




The integration of Eq. (52), the details of which are 
presented in Appendix E, yields, 




(53) 



for t < 0, and t > k 



and 




2 



(5*0 



for 0 < t < k 



The deflection angle B can be found by letting t in 

Eq. (53) j or by letting t in Eq. (50) and then performing 

the integration as shown in Appendix F. The deflection angle 
B is found to be uniquely related to the parameter k as 



To obtain the shape of the bucket and -the actual location of 
the free streamlines, Eqs . (53) and (5*0 are replaced in 
Eq . (*I5) and integrated. 

The integration must be performed numerically using the 
method described by Hartree [20] because of the singularity 
at the point t = 1. It must be pointed out that the 
singularity at t = 1 is an integrable singularity of the form 



B V + 2 (k - 1) 



[■ 




^Tk 



-1 (55) 
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1 



dx 



S 

0 v/T 

No attempt is made here to perform the numerical integration 
of Eq. (^5) even though the procedure is well known as shown 
in one application by Tinney [2], because the method previously 
described, namely, Levi-Civita ' s method, is very similar to 
the one presented herein and has the advantage of being very 
flexible to obtain different bucket shapes and presents no 
singularities for the integration. 

D. A CRITIQUE OF THE ANALYTICAL METHODS FOR 

TWO-DIMENSIONAL FLOWS 

The past as well as the present analyses have emphatically 
demonstrated that almost all of the difficulties of the 
analysis of jet deflection from curved obstacles are ascribable 
to the determination of the initially unknown free streamlines. 
In view of this fact a direct analysis of the problem, that is, 
the determination of the jet deflection and the position of 
the free streamlines via either purely analytical or partly 
analytical and partly numerical techniques is at best an 
exceedingly complex proposition. As propounded by Birckhoff 
et al . [ 1 4 3 and other hydrodynamicists what one can hope for 
is an indirect solution of the problem; that is, obtaining a 
family of obstacle shapes which will yield the prescribed jet 
deflection angle and provide a familiarity between the shapes 
so obtained and the range and change of direction of the free 
parameters involved in the analysis. What has been said so 
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far holds true for both the plane and axisymmetric flows. 
Should one adopt a purely numerical procedure (finite differ- 
ence, relaxation, finite element, etc.) it is then and only 
then that it is possible to begin with a prescribed obstacle 
shape and to approach a unique solution through successive 
iterations . 

In what follows, attention will be focused only upon the 
evaluation of the individual as well as relative merits of 
the analytical methods dealing with the deflection of plane 
jets . 

The hodograph-method developed herein is capable, within 
the limitations of the potential-flow theory, of analysing 
the flow in curved obstacles and generating a family of 
practically suitable bucket shapes. Furthermore, this method 
has the inherent capability of being extended to the genera- 
tion of even a larger family of obstacle-shapes through the 
use of distributed sources and sinks in the hodograph-plane . 
Even though this concept has not been pursued further herein, 
there are no conceptual difficulties in doing so. 

The method devised by Levi-Civita and successfully em- 
ployed herein requires the judicious selection of the varia- 
tion of 6 with 0 and the use of a suitable series to enable 
one to perform the integrals leading to the coordinates of 
the obstacle as well as of the free surfaces . The deflection 
angle may be determined either by calculating the asymptotic 
slope of the free streamlines or by finding the pressure 
distribution over the bucket and using the equation of 
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momentum. The emerging fact is that a brief familiarity with 
this method enables one not only to generate a family of 
suitable curved obstacles but also to round the corners of 
buckets composed of straight segments. For example, by 
letting (a^ - a * 0 one finds that the bucket shape 
generated by this method exactly approaches that previously 
solved by Sarpkaya [1] (see Figure. 8a). For slightly 
larger values of (a^ - a 2 ) the corner tree Fig. 7) is 
rounded, thus enabling one to evaluate the effect of the 
deviation of the segment PQ from a sharp corner on the 
deflection angle. In fact, the results show that for a 
given b/a and a/d ratios (see Fig. 1) the effect of the 
curvature is relatively small as anticipated by Sarpkaya [1] 
on the grounds that the flow of a real fluid near the sharp 
corners form nearly stagnant regions comprised of small 
corner vortices and provides a slightly rounded streamline. 

Finally, it should be noted that Levi-Civita ' s method 
has been rendered practicable through the introduction of 
the Fourier-series into the analysis. It may be easily 
extended to the analysis of the cases where the lip angle of 
the bucket is other than tt or the bucket has a cusp on its 
axis of symmetry or to the cases where the impinging jet is 
not free and emerges from a prescribed nozzle. Such shapes 
are known to be used in the design of spillways, Pelton 
wheels, two-dimensional cascading thrust-re versers , etc. 

As to the Riemann-Hilbert method, it may have already 
been evident that it is in essence a mathematically more 
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elegant and sophisticated version of Levi-Civita ' s method in 
prescribing certain functions. It is for this reason that 
no particular solutions have been presented and only a linear 
variation of f(t) was studied. Needless to say, a quadratic 
or cubic variation of f(t) could have been employed. 

Evidently the resulting integrals would not have been 
integrable in closed form. Finally, it should be pointed 
out that the Rlemarm-rTilbert method is perhaps the only 
method aside from purely numerical methods with which one 
can incorporate into the analysis the effect of body forces 
such as gravity. 

In closing the discussion of the methods dealing with the 
deflection of plane jets from curved obstacles, one must once 
more emphasize that all of the foregoing analyses and solu- 
tions deal with plane Laplace-field problems where the real- 
fluid effects must be out of necessity ignored. Because of 
the limitations imposed upon man's understanding of turbulence, 
the development of laminar and/or tui’bulent boundary layers 
under a nonzero pressure gradient along curved walls, entrain- 
ment between the jet and the surrounding fluid which may or 
may not have the same physical properties, and finally the 
Coanda effect (the tendency of jets to attach to adjacent 
surfaces because of the entrainment deprivation) , the real- 
fluid effects on jet deflection cannot as yet be analyzed. 

It is because of these reasons that the analysis of the 
problem with the assumption of an inviscid fluid becomes 
more important for It Is only through the comparison of such 
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solutions with the experimantal results that one can assess 
the degree of significance of the real fluid effects on the 
governing parameters. As stated in the introduction and as 
verified by numerous experiments (Robertson [13] , p. 626) > 
the phenomenon of jet deflection is governed primarily by 
pressure and inertia forces and that the real-fluid effects 
are of minor significance in determining the forces acting 
on the body and the deflection of the jet. . These facts will 
be stressed again in connection with the discussion of the 
analytical and experimental results for axisymmetric jet 
deflection to be subsequently presented. 
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III. ANALYSIS OP AXISYMMETRIC JET IMPINGEMENT 



A. INTRODUCTION 

As stated in the general introduction, the two- 
dimensional counterpart of the jet-deflection problem has 
been treated by several investigators through the use of 
numerical and approximate analytical methods. Additional 
techniques based on the modified-hodograph method, Riemann- 
Hilbert integral transforms, and the Levi-Civita* s method 
have been successfully developed and employed in the 
preceding section in the analysis of the jet impingement 
upon two-dimensional curved buckets. 

The three-dimensional counterpart of the jet-deflection 
problem has not yet been solved in any generality. Attempts 
to formulate an exact solution have been mostly unsuccessful 
even for the axisymmetric inviscid flows with no body forces. 

The case of a circular jet striking a plate normally was 
analyzed by Schaach [21] using approximate methods similar 
to those of Trefftz [22] with successive adjustment of the 
free streamlines and through the use of the Fredholm integral 
equations. These results are in good agreement with those 
obtained experimentally and to the surprise of the many 
investigators, they are almost identical to those of the 
two-dimensional cases. 

Other noteworthy contributions to the analysis of the 
jet-efflux from nozzles and orifices with straight boundaries 
were made by Southwell and Vaisey [23], Rouse and Abul-Fetouh 
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[ 2 * 1 ], Garabedian [25], and Hunt [ 26 ] through the use of the 
relaxation and finite difference techniques, and singular- 
integral equations. Schnurr et al . [27] used the relaxation 
method to analyze the turning of two-dimensional and axisym- 
metric jets from curved surfaces where there was only one 
free stream surface; i.e. the jet was assumed to leave the 
deflector exactly parallel to the tangent at the lip of the 
deflector surface with a uniform velocity distribution and 
the consequences of the difference between the actual deflec- 
tion angle and the said tangent to the deflector surface was 
taken as a measure of the "spillage" and expressed in terms 
of a "turning-effectiveness coefficient" determined experi- 
mentally. Thus, this analysis does not constitute a solution 
to the problem under consideration. 

Jeppson [ 28 ] proposed an inverse formulation of the 
problem for the flow exiting from a circular orifice by using 
x and r as dependent variables and <{> and \p as the independent 
variables. Through the use of this technique Jeppson [29] 
also analyzed the flow of a jet from a convergent nozzle with 
curved boundaries. 

All of the investigations described so far have attempted 
to provide better solutions to the various configurations of 
the orifice flow. These techniques suffer in general from 
convergence and accuracy problems and nearly all resort to 
simple trial-and-error procedures to locate the free surfaces 
and to satisfy the boundary condition that the free surfaces 
be streamlines of constant velocity. Suffice it to say that 
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a method was needed which could yield solutions of a 
prescribed accuracy for a wide variety of fluid flow problems 
involving Dirichlet, Neumann, and mixed boundary conditions. 

Zienkiewicz and Cheung [30] proposed in 1965 the applica- 
tion of the finite element method to the solution of field 
problems involving the equations of Laplace and Poisson. 

Since then a significant number of applications of the finite 
element method to fluid dynamics has aooearsd in the litera- 
ture (see e.g. Norrie and de Vries [31] and the references 
cited therein) . 

The method has been recently applied to several jet-efflux 
problems involving only one freestream surface and relatively 
small jet contraction by Chan and Larock [32]. Suffice it to 
say that the finite element method has proven its versatility 
and applicability to a wide variety of two-dimensional and 
axisymmetric f low-problems . The truly three-dimensional 
problem such as the oblique impact of a circular jet on a 
plane surface, with or without gravitational forces, just to 
name one, offers considerable challenges and demands additional 
sophistication. 

The special case of the gravitational force acting along 
the axis of symmetry may easily be analyzed by modifying the 
boundary conditions such that the velocity along the free 
surface varies according to the Bernoulli theorem. 

This part of the investigation encouraged by the advent 
and the continued success of the finite element method is 
devoted to a determination of the angle of deflection, the 
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location of the freestream surfaces and the velocity and 
pressure distributions caused by a finite curved boundary 
in general and a hemispherical boundary in particular placed 
symmetrically with respect to the axis of an axisymmetric, 
inviscid jet issuing from a nozzle. This problem has not 
been previously analyzed through either exact or approximate 
methods . 

As stated earlier, the deflection o£ a jet by a solid 
boundary is well suited to potential-flow analysis because 
of the dominance of inertia and pressure forces in the 
establishment of the flow pattern. The result in general 
represents asymptotic values which are approached as the 
effect of secondary variables such as entrainment, boundary 
layer, compressibility, jet attachment to adjacent surfaces 
(Coanda effect), etc. are decreased. Specifically, if the 
total angle through which the jet is deflected is determined 
for idealized conditions, then the principle of impulse and 
momentum can be used to compute pressures, forces, or other 
dynamic characteristics of the system such as the reverse- 
thrust ratio. In this manner one can not only optimize the 
geometry of the system for optimum performance, measured in 
terms of reversed thrust, but also can assess the significance 
of the real-fluid effects, cited above, on its performance. 

B. ANALYSIS 

The problem may be defined in general with reference to 
Figure 10. An axisymmetric jet issues from a nozzle BC 
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Figure 10 General case of axisymmetric bucket an nozzle 



with an entrance velocity V Q at AB. The jet impinges symmet- 
rically upon the curved bucket LKH and then the deflected jet 
leaves at an angle 3- The upper and lower surfaces of the 
jet; namely, CDE and HG, asymptotically converge to the 
centerline of the jet in accordance with the conservation of 
mass . 

It is apparent from Fig. 10 that for a given nozzle and 
bucket shape, the deflection angle varies with the diameters 
of the bucket and the nozzle and the separation distance 
between the lips of the nozzle and the bucket. In other 
words, 3 can be expressed as 
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(56) 



3 = 3(R 0 /r 



o 



s/r ) 
o 



Thus, the determination of this functional relationship 
constitutes the basis of the problem. 

The fluid is assumed to be inviscid, incompressible, and 
free from body forces. Consequently, the problem may be 
formulated in terms of either the velocity potential <{> or 
the Stokes stream function ip. Here the potential function 
is chosen as the unknown field parameter. 

The governing equation is then given by 



ifi . i m + 

ax 2 r 




(57) 



in which x and r are the axial and radial coordinates. 

Because of the axial symmetry there is no variation in the 
meridional direction. This equation is to be solved with 
the boundary conditions which simply states that the normal 
component of velocity be zero along both the solid and free 
surfaces, i.e. the free surfaces be streamlines of constant 
velocity . 

For an axisymmetric flow, the solution to the Laplace- 
field equation satisfying the specified normal-velocity 
boundary condition (3<J)/3n) a is given by the admissible 
function <J> which minimizes the functional (see e.g. Zienkiewicz 
and Cheung [30]) 

I(<t>) = pit // [(f^-) 2 + (f^O 2 ]rdxdr - 2px 6 <t>(|4) a rdS (58) 

^ d X dr q d n 
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in which A is the half of a meridional section of the flow 



and C is a portion of the curve bounding this area where the 
normal derivatives are prescribed. The first integral in 
Eq. (58) represents the kinetic energy of the fluid within 
the entire control volume and the second integral represents 
twice the work done by the impulsive pressure p4> on the 
boundaries in starting the motion from rest (Milne-Thompson [ 4 ] ) . 

It is apparent from the foregoing that in the finite- 
element method the boundary conditions are an integral part 
of the functional 1(0) and will not have to be dealt with 
separately as in the case of finite-difference or relaxation 
techniques. This is, in fact, where the beauty and the power 
of the finite element technique lie. 

The approximate minimization of the functional 1(0) is 
accomplished as follows. Firstly, the field of interest is 
divided into M triangular elements and a functional 1(0) is 
written as a sum of M element-functionals I (0) each of which 
has the form of Eq. ( 58 ), i.e. 



M e 

1(0) = E I e (0) 
1 



(59) 



Figure 11 shows the numbering system of the corners and mid- 

0 

points and the corresponding area coordinates = A^/A 
in which A^ is the area of a subtriangle and A e is the area 
of the entire triangle. Then it is clear that + ^2 + ^3 = '*'‘ 
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3 (0,0, I ) 




Figure I! A triangular element with area coordinates 



Secondly, a polynomial of sufficient order is selected to 

0 

represent the variation of the potential function <j> within 
each element. 

Apparently, a linear variation of <j> e leads to constant 
velocities along a given line within an element and to 
unacceptably large errors in the velocities between neighboring 
elements. A second-order polynomial, however, allows at least 
a linear variation of velocity components. 

The use of higher order polynomials may increase the 
accuracy of the results and the convergence of the minimization, 
but it may also lead to extremely tedious arithmetical manipu- 
lations. For these reasons and following Felippa [33] and 
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Chan and Larock [32], a second-order polynomial has been 
chosen. In terms of the area coordinates, this leads to 



< !> e = <f>J Xi (i = 1,. . . , 6) (60) 

0 

where <j>^ are the coefficients of the polynomial and corresponds 

0 

to the values of <j> at the node i. It is important to note 
0 

that <j>^ must be regarded as constants during the integration 
of Eq. (58) in each element. 

The term y^ may be written after Felippa [33] as 



x i = CC 1 (2^ 1 -1) ,C 2 (2^ 2 -1) ,C 3 (2c 3 "l) } ^C 1 C2^C 2 C 3 > J k 1 C 3 ] 
The velocity components u and v are given by 



u 



_ 

9x 



and 



v 



31 

9r 



( 61 ) 



Then the velocity in each element can be written as 




3Xi 

9x 



(i = 1 to 6) 



( 62 ) 



v e = <j>? - 5 —= - (i = 1 to 6 ) ( 63 ) 

Y i 9r 

The evaluation of Eqs . (62) and ( 63 ) require the derivatives 
of y^ with respect to x and r or the evaluation of terms 
such as 
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3 5- 



i 1 = (4c 1 -l) — = (^5,-1) -— = T® 

3x - 1 3x x 2A X 



r 5 ? (25 -l)n 3 5„ b„ 
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2A 
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in which 



3q 

ax 




and 



b 



k 



r 



0 



The foregoing closely follows the notations developed by 
Felippa [33]. 
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It is evident that the derivatives of x with respect to 
r may easily be evaluated by changing x to r or b^ to a^ 
given by 

a k = x j " X i 

and denoting the resulting array of 1\ by 1\ . Thus, one 
can write 




(4^. - l)b . 

s i a 



(no sum on i) 



(65a) 




(65b) 



and likewise 




(no sum on i) 



(66 a) 




( 66b ) 



In which 



2A 



e 




and 



r = + r 2 ; 2 + 
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The minimization of the functional I ( <f> ) in accordance with 
the Ritz method requires that the partial derivatives of 



I ( 4> ) with respect to be zero. In other words 



ai e (40 = o 



(67) 



94>. 

Y i 

which may be evaluated in a manner quite similar to that 
used in Eq. (6^1). Thus, Eq. (67) reduces to 



9l e (4>) 

94>? 



2 pir // (T?T® +T?T!)<}>?rdA - 2irp (fi x, (l^) a rdS = 0 (68) 

ljijj pi on 

A e L 



or to 



a J e ( U . _ K ? 

3<f e ^ •> 

Y 1 




(69) 



in which and B*T represent, respectively, in the structural- 

mechanics terminology, the element stiffness matrix and the 
load matrix for a triangular element. In fluid mechanics 
terminology, they simply represent the geometrical configura- 
tion of the grid within the domain and the boundary conditions 
imposed on it. These coefficients have been previously 
obtained by Chan and Larock [32] and are tabulated in Appendix 
G , 

In fluid-flow problems where one or more parts of the 
boundary are to be determined as part of the numerical 
analysis, the solution of the system of equations represented 
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by Eq. ( 69 ) through the use of the Gaussian elimination 
procedure requires the use of a suitable (i.e. accurate and 
convergent) iteration scheme. This important phase of the 
analysis will be described in the following. 

C. ITERATION SCHEME 

The iteration scheme consists of the designation of a 
moving and fixed grid, the adoption of a scheme for assuming 
the initial position of the deflected jet, and the successive 
adjustment of the free surfaces until the boundary conditions 
are satisfied. 

1 . Designation of a Grid 

Evidently, the minimization process cannot be carried 
out with a fixed grid as it would be in a structural system 
subjected to small deformations, primarily because the 
initial position of the free surface is not a priori known, 
and the initially-assumed boundaries may be too far off from 
their final positions. In view of the foregoing, the fluid 
domain was divided into a moving and a fixed grid (see Fig. 
12). The fixed grid occupies that region of the flow which 
is not, in the final analysis, to be intersected by the free 
surfaces. This does not impose any restrictions on the anal- 
ysis and requires no more than a passing familiarity with 
the elementary principles of jet deflection from solid 
boundaries. Should the boundary of the fixed grid be inter- 
sected by the free surface during the iteration process, a 
new fixed grid has to be assigned. 
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Figure 12 Jet deflection from a hemispherical bucket (R /r =1.6, s/r -0.8) 



The moving grid is comprised of triangular elements 
whose nodal points may be moved along prescribed lines prior 
to the commencement of the next iteration. This procedure may 
best be described with reference to Pig. 13. At the end of 
a given iteration, the nodal points A, B, C are moved inward 
or outward, in a manner and magnitude to be described later 
and the positions of the nodal points F, G, E, H, D, and K 
are recalculated such that 

AF = 0.15 AN 



and 



AG = 0.50 AN 

and similarly for BE, BH, etc. Evidently, the coefficients 
0.15 and 0.50 are arbitrary and different values could have 
been used provided that the resulting elements maintain a 
shape more or less compatible to their original shape and 
that the ratio of the sides of a given element does not 
significantly differ from unity. 

The points intermediate to AB and BC, i.e. the points 
1 through 6, are calculated at each iteration through the use 
of a subroutine called CURVE-FIT. This subroutine is a third 
order polynomial connecting the points A, B and C, and evenly 
spacing the points 1 through 3 between AB and the points 4 
through 6 between BC. It should be immediately pointed out 
that there is no need for an exotic subroutine to perform 
this procedure and that any other judiciously-selected 
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(b) Moving grid in the deflected jet 



Figure 13 Construction of the moving grid 
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polynomial could have been employed. The sole purpose of 
the foregoing procedure is to smooth the boundary between 
two successive iterations and to obtain a better representa- 
tion of the curved streamlines. Furthermore, this procedure 
eliminates the occurrence of localized low or high velocities 
which could have stemmed from a zig-zag positioning of the 
points such as 1, 2, 3 and thereby accelerates the convergence 
process. In the deflected portion of the free jet (see Fig. 
13b) the points such as P and R are moved inward or outward 
as dictated by the iteration conditions and the point Q is 
located in the middle of PR prior to each new iteration. 

2 . Assumption of Initial Boundaries 

The assumption of an initial boundary requires the 
consideration of the equation of continuity, some familiarity 
with jet deflection problems in fluid mechanics, and possibly 
the use of a deflection-angle predictor, if one is available. 

Ordinarily, a deflection angle 3 and the length of 
the deflected free jet (HG) are assumed (see Fig. 14) . Then, 
through the use of V , r , r„, and the equation of continuity, 

O O r 

the thickness d of the jet is calcuated, to a first order of 
approximation by 



d 



V r 
o o 



2r 



F 



(70) 



Subsequently, the points G and E are located such that F is 
intermediate to both. Then the point G is connected to H 
(the upper lip of the jet) by a straight line, and the point 
E is smoothly joined to the lower portion of the assumed 




Figure 14 Nomenclature for a hemispherical bucket and 
straight-nozzle combination 



boundary. In passing, it should be noted that the approximate 
use of the equation of continuity in sections intermediate to 
GH and ED may significantly simplify the assignment of the 
initial position of ED. For this purpose, all one has to do 
is to assume a uniform velocity at each section and calculate 
the thickness of the jet at the section under consideration 
through the use of an expression similar to that given by 
Eq. (70). 

The description of the assignment of the fluid domain 
will not be complete without a brief discussion of the length 
of the nozzle and the deflected jet. Clearly, the sections 
AB and EG must be so located that one can for all intents and 
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purposes assume the velocities at the said sections to be 

uniform. Otherwise, the analysis will not be correct and the 

jet boundary will be influenced by the conditions created 

by the proximity of either the section AB and/or EG to the 

points C and H. To this end, preliminary calculations were 

carried out with nozzle lengths (BC) as large as 6r Q and it 

was found that it could easily be reduced to 3r Q without 

materially affecting the accuracy of the velocity at AB. In 

fact, the calculated velocities at AB with a nozzle length 

of 3r were accurate within ± 0 . 5% * Calculations were per- 
o 

formed with a nozzle length of 3.8r o ,i.e. about twice the 
bucket diameter. 

The deflected- j et length GH was varied between 0.5r o 
and 0.7r o depending more or less on the particular angle of 
deflection. Suffice it to say that no noticeable change was 
observed in the calculated position of the jet. This is 
primarily because of the fact that for a jet length of 0.6r o , 
the ratio of the jet length to final jet thickness is about 
3.0 and that according to all the previous calculations and 
exact solutions performed for two- and three-dimensional jets, 
the jet attains its asymptotic uniform velocity profile in 
a distance less than 3d. 

3 . Iteration of the Boundaries 

In the foregoing, all of the necessary elements of 
the analysis and the domain have been provided for the com- 
mencement of the iteration process. We will now describe 
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the process which eventually leads to the final position 

of the jet for a given geometrical positioning of a 

hemispherical bucket and a straight nozzle (see Fig. 14) . 

To begin with, the free surfaces were assumed to 

be solid boundaries. Thus, the two boundary conditions; 

namely, constant pressure and zero normal velocity along the 

free surfaces, reduces to a single condition. Specifically, 

the freestream surfaces ^should be streamlines of constant 

velocity. The normal velocity was written zero everywhere 

along all the rigid and free surfaces except along EG (see 

Fig. 1*1). An arbitrary value of <t> = 100 was assigned along 

EG. Needless to say, the final solution does not depend on 

this particular value of <j> . 

To start a computation the velocity V q along AB was 

assumed to be equal to unity and the velocity V was calculated. 

c 

For an exact solution V should be equal to unity assuming 

c 

that the velocities along all the free surfaces are taken 
equal to unity. Evidently, V will turn out to be larger 

G 

than unity because of the fact that the pressure along AB is 
larger than zero and thus the fluid accelerates along BC and 
decelerates along AL. The velocity at C was immediately 
corrected to a first order of approximation by iterating on 
V q through the use of the following FORTRAN STATEMENT: 

V0 = V0 - (VC - 1.0) * V0 (71) 
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This process of iteration has been used in all subsequent 
iterations. Apparently, the iteration on could have been 
accelerated by beginning the analysis with a V Q value smaller 
than unity (say, V q = 0.9)* 

Subsequently, the velocities along the free boundaries 
were calculated everywhere through the use of the <f> values 
at the nodal points on the free surface and, as anticipated, 
it was found that they were not everywhere equal to unity. 

The reduction of these velocities to unity within a prescribed 
error constitutes the essence of the iteration process and 
the use of the moving grid system. 

Let V s represent the velocity at a point along the 
free surfaces. Then, the free surface at that point was 
moved inward or outward, depending on whether V was greater 

3 

or smaller than unity, by an amount given by the following 
statement : 

DELV = (VS ** 3 - 1.0) * FAC (72) 

in which FAC is a multiplier. The correction DELV was applied 
to the nodal points along a line normal or nearly normal to 
the boundary. Various values of FAC have been used in the 
preliminary calculations, and it was found that a FAC value 
smaller than 0.05 will yield a convergent iteration. In the 
present analysis either one of the following two techniques 
have been used in correcting the free surface. In the first 
method a single FAC multiplier of 0.015 was applied to both 
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the upper and the lower boundary, five iterations were per- 
formed, and then the angle of deflection 3 was increased or 
decreased by two degrees depending on whether the velocities 
on the upper free surface was smaller or larger than unity. 
This process has been continued until everywhere the veloci- 
ties reached a value of 1 ± 0.01. It is important to note 
that the angle of deflection, the jet thickness, and the 
radial coordinate end r^ change during each iteration pri- 
marily because of the corrections made on the free surfaces. 
However, the change in 0 per iteration was about 0.25 degrees 
Consequently, one would have needed at least 40 iterations to 
reach a correct jet position from an initially assumed posi- 
tion with an error of ± 10°. Evidently, then, the ± 2° 
increments applied to 3 have accelerated the convergence 
process . 

In the second method two FAC values were used; namely 
FACU = 0.001 and FACL = 0.015. Five iterations were carried 
out by applying FACU to the upper portion of the jet (along 
GH in Fig. l J t) and FACL was applied to the lower boundary. 

At the end of the fifth iteration the role of the FAC values 
were interchanged and five more iterations were carried out. 
This simply amounted to first iterating on the lower boundary 
and then on the upper boundary. As before, the deflection 
angle was incremented by two degrees until the correct jet 
deflection angle has been arrived at, as evidenced by the 
uniformity of the velocities along the upper and lower free 
surfaces. 
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The full description of the iteration process would 
not be complete without stressing few other computational 
procedures which have enabled the prediction of the initial 
deflection angle and signaled at the arrival of the correct 
jet position. These will now be described. 

The boundaries were iterated upon through the use 
of either one of the two methods cited above until a deflec- 
tion angle has been arrived at which was not more than a few 
degrees off from the final jet position. At that time no 
limitations were imposed on the number of iterations and the 
jet was allowed to arrive at its final position. In passing, 
it should be noted that the programs were written in FORTRAN 
IV and double precision arithmetic was used. Since the pro- 
grams on the Naval Postgraduate School computer (IBM 36 0 — 67 ) 
had to be run overnight for a period of about one hour, the 
stops between two successive runs have provided the writer 
with an opportunity to examine the variation of the veloci- 
ties, to observe the speed of convergence and finally to 
decide as to whether one has arrived at the stage of letting 
the computer finally carry out twenty to twenty-five itera- 
tions to establish the final jet position. 

Ordinarily, in dealing with the analysis of jet de- 
flection problems from arbitrarily-shaped curved surfaces, 
one will have at first very little insight as to the value 
of the deflection angle which the jet will finally assume. 
However, after one or two calculations with one or two geo- 
metrical combinations of a given bucket and nozzle shape. 
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one will not only be able to guess fairly accurately the 
initial jet-deflection angle but will also be able to derive 
a predictor equation. We will now describe the derivation 
of such an approximate equation which has been useful in 
predicting an initial deflection angle for all the remaining 
cases following the study of the first case. 

Consider the equation of momentum applied to the 
fluid domain under consideration: 



P nr + Trr pV 
o o o K o 



t j.- n v v . 1 
O K O J 



in which P represents the uniform pressure along AB, V. the 
o j 

velocity along the free streamlines (V. = 1 in the present 
analysis), p the density of fluid, and T the force required 
to hold the bucket at rest. Evidently, if V^, (3, and 
were known, then the force acting on the bucket could have 
been calculated. In fact, the present analysis does enable 
one to do so not only through the use of the momentum equa- 
tion but also through the integration of P dA (dA is the 
projected elemental area of the bucket at the center of which 
the pressure is P). Since the velocities along the bucket 
were known for the first case, P dA was integrated and the 
resulting thrust was calculated. This procedure has not 
only provided a means of checking the thrust calculated from 
the momentum equation, but also an equation to predict the 
deflection angle for other geometrical combinations of the 
nozzle and the bucket. Thus, writing 
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/ PdA = 2 r 



o 




j j 



one has 



r 



~2 / PdA - tt[1 + (V 0 /V j ) 2 ] 



cos 3 



o 



(75) 



2tt(V o /V.) 



The equations so far derived are exact and Eq. (75) could be 

used as a mean to predict the angle of deflection for other 

R /r and s/r values if one assumes that / PdA remains 
o o o 

essentially constant. Thus using the / PdA value obtained 
from the first case and a V q value of approximately 0.8, 
one arrives at the initial value of 3- It should be stressed 
that the foregoing procedure is designed simply to help the 
analyst to arrive faster at the correct solution. It neither 
alters the final result nor is absolutely necessary for the 
solution of the problems discussed here. 

D. HEMISPHERICAL TARGET-TYPE THRUST REVERSERS 

.. The use of thrust reversers on modern jet craft for the 
purposes of decreasing the landing load, maintaining 100$ 
engine speed vjhile landing, decreasing the thrust while 
maneuvering, reducing the length of the landing strip, tire 
wear and tear, etc. has been one of the most important 
developments in aircraft industry in the last two decades. 

To perform effectively, a thrust reverser should yield a 
relatively small jet-deflection angle, should have easy 
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stowage configuration, should have relatively small boattail 
or stowage drag, and finally should not adversely affect the 
engine performance. A hemispherical thrust reverser appears 
to fulfill these requirements. Furthermore, there exist 
some experimental data on this type of reverser (Steffen et 
al. [3^] and Povolny et al. [35]) for comparison with the 
results of the present analysis. 

So far as the finite element method is concerned, there 
is no restriction on the nozzle and the bucket shape and they 
could be of any axisymme trie , reasonably-chosen configurations. 
In other words, the finite-element method enables one to 
optimize, within the assumptions of the inviscid flow theory, 
the relative position of a given bucket and nozzle, the effect 
of some geometrical modifications on the bucket and/or nozzle 
shape (such as placing a flat plate within the bucket) on the 
performance of the reverser. Furthermore, it enables one to 
assess the degree of significance of the real-fluid effects 
such as viscous and turbulent energy dissipation, entrainment, 
Coanda effect, nozzle-pressure ratio, etc. on the actual 
performance of the bucket. 

The geometrical configurations, namely, the ratios R 0 / r 0 
and s/r Q were chosen partly on the basis of the past experi- 
ence with such reversers. For example, it has been experi- 
mentally demonstrated by Steffen et al . [3^] that the larger 
bucket, the larger is the reversed thrust and that the closer 
the nozzle to the bucket, the larger is the back-pressuring 
effect. In other words, when a nozzle is placed very close 
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to the bucket, the pressure rises in the nozzle because of 
the blockage of the flow by the bucket and this, in turn, 
tends to affect the engine performance. Therefore, there is 
an optimum nozzle-bucket configuration which yields optimum 
reverser performance. 

It is on the basis of such considerations that the nozzle- 
bucket geometries listed in Table I have been chosen for the 
numerical analysis and experiments. 



TABLE I 




Geometrical Charact 
of the Hemispheric a 
buckets analyzed 


1 ..... 


Angle of Deflection 3 (deg. ) 


L 

Finite 


element analysis 


Experimental 

(water) 


Case 

No. 


R /r 
o o 


slope of 

the inner 

s/r surface 

o 


mean 

deflection 

angle 


slope of 
the outer 
surface 


slope of 
the outer 
surface 


1 


1.8 


0.4 


21.1 


20.3 


19.5 


18 


2 


1.8 


0.8 


37-0 


35.0 


33.0 


30 


3 


1.6 


'0.8 


53.2 


50.5 


47.8 


46 



The coordinates of the free surfaces for the three cases 
cited in Table I are presented in Appendix H. The computer 
program used in the study of the case number 3 is given in 
Appendix I. Suffice it to say that similar programs were used 
for the analysis of the other two cases. The miner differences 
between the three computer programs essentially resulted from 
the use of different moving and fixed grids to suit the 



particular geometry under study. An effort has been made to 
keep the number of elements and the bandwidth in the resulting 
matrix as low as possible and yet maintain a reasonable 
degree of accuracy. At various stages in the evolution of 
the computer programs larger or smaller number of elements 
have been used, and it was found that the number of elements 
used herein (288) is quite sufficient. Needless to say, the 
coordinates of the fr’se surfaces may be calculated to a 
greater degree of precision by increasing the number of the 
elements particularly in the moving grid. The additional 
corrections to be applied to the nodal points on the free 
surfaces at the end of the final iteration were less than 
10 . It is estimated that the calculated deflection angle 

will not vary more than il.0°. In view of the fact that the 
real-fluid effects such as entrainment, Coanda, etc., can 
decrease the deflection angle by as much as ten degrees, it 
was not deemed necessary to pursue the calculation of the 
deflection angle to a greater degree of precision. 

A series of experiments were conducted with water through 
the use of straight, sharp-edged nozzles of appropriate 
dimensions (see Table I) with a hemispherical cup of 0.72 
inch radius primarily for the purpose of determining the 
slope of the outer surface of the deflected jet. The slope 
was determined from photographs similar to that shown in Fig. 
15 and tabulated in column 7 of Table I. 

It is evident from Table I that the experimentally deter- 
mined values of 3 are somewhat smaller than those obtained 
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numerically (see columns 6 and 7, Table I). The reasons for 
these differences are several, the most notable ones being 
the energy dissipation due to viscous action and the Coanda 
effect. The latter will be discussed in detail further later 
in connection with the experiments conducted with air. 

Suffice it to say that the Coanda effect with water is 
relatively small and the difference between the computed and 
the measured angles of deflection is within the range of 
experimental errors and numerical approximations. 

The theoretical and experimental results may also be 
compared in terms of a measure of the reversed- thrust because 
of its obvious practical significance. 

In the aircraft industry the performance of the reverser 
is expressed in terms of a reverse-thrust-ratio defined by 

actual reversed thrust 

forward thrust of the nozzle alone 



or in other words by 



n R = 



it r p V. V cos £ 
o .1 o 

? 2 
7 T r ^ p V . 
o H 0 



V 

O n 

— COS 3 



The n D values found in the present analysis, and those 
n 

obtained experimentally by Steffen et al. [34] are tabulated 
below . 
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TABLE II 



Case 

No. 


n R 

Present Analysis 


n R 

Experimental [34] 


1 


0.766 


0.82 


2 


0.725 


0.79 


3 


0.515 


0.62 



Evidently the ri R values obtained experimentally by Steffen et 
al. [34] are somewhat larger than those obtained in the present 
analysis. As noted earlier, the observed differences are 
primarily due to the Coanda effect and the nozzle-pressure 
ratio. The Coanda effect decreases the deflection angle and 
thereby increases the reversed thrust. In the experiments 
conducted by Steffen et al . [34] the outer surface of the 
nozzle was streamlined in the form of a boattail to decrease 
entrainment and thus to decrease 8 or to increase In 

fact, as noted by Steffen et al . [34]» the pressure reduc- 
tions on the boattail were large enough to account for as 
much as 10 to 15 degrees reduction in the deflection angle. 

Thus the relatively large differences between the values 
computed in the present analysis and those obtained by Steffen 
et al. [34] are primarily attributable to the Coanda effect. 
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In order to illustrate the influence of the Coanda effect, 
one particular experiment was carried out with air using a 
1.25 inch radius nozzle and a 2 inch radius hemispherical 
bucket with a 1 inch spacing (i.e. with R o /r Q = 1.6 and 
s/r Q = 0.8). The photograph of the deflected jet (visualized 
with aluminum pov:der) is shown in Fig. l6a. Figure 16b shows 
the deflection of a water jet from a similar nozzle-bucket 
geometry. It is evident from the comparison of these two 
photographs that whereas the slope of the outer free surface 
with air is 33 degrees, with water it is H6 degrees. The 
13 degree reduction in the deflection angle due to the Coanda 
effect corresponds to an approximately 20$ reduction in . 
This brings the calculated rip values into closer agreement 
with those experimentally obtained. 

The nozzle-pressure ratio (total nozzle pressure/ambient 
pressure) or the actual velocity in the nozzle causes varia- 
tions in 3 primarily because the entrainment needs of the 
deflected jet and hence the Coanda effect Increases with 
increasing jet velocity. 

Suffice it to say that the analysis of thrust reversers 
is extremely complex for the cases where the jet impinges 
asymmetrically or where the gravitational effects are impor- 
tant. Furthermore, some of the practical problems associated 
with thrust reversal, such as the reattachment of the jet to 
the nacelle of the engine, hot gas reingestion, interaction 
of the deflected jet with the ambient stream and the adjacent 
surfaces (in this case the outer surface of the nozzle) cannot 
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(b) 



Water jet 




Figure 16 



Air and water jets deflected by a 

bucket (R /r = 1.6, s/r = 0.8) 
o o o 



hemi spheri ca I 
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yet be analyzed by any one of the existing numerical methods, 
even for relatively simple two-dimensional geometries 
composed of straight boundaries. Entrainment depends on the 
velocity distribution and the turbulence level in the jet, 
shape of the adjacent surfaces, the physical properties of 
the jet and of the fluid medium into which the jet discharges 
Most of these parameters are obviously interdependent and at 
present there does not seem to be any hope that their effect 
on the deflection of an axisymmetric jet can be evaluated. 

Nevertheless, approximate analyses such as the one 
presented herein help to isolate the more promising types of 
thrust reversers, investigate their potentialities, evaluate 
their ideal performance characteristics, and to delineate 
the range of importance of the geometrical variables involved 
Then those cases that are shown to be promising through 
computer experiments can be tested in a laboratory, at a 
sufficiently reduced expense, to evaluate the effect of 
additional fluid-mechanical phenomena on their actual 
performance . 
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IV. SUMMARY OP RESULTS AND CONCLUSIONS 



A class of analytical and numerical solutions to the 
problem of inviscid-jet deflection from plane and axisym- 
metric concave surfaces have been presented. The methods 
of solution for the two-dimensional cases included the 
modif ied-hodograph method developed herein as well as those 
devised by Levi-Civita and Riemann and Hilbert . 

The modified-hodograph method is capable, within the 
limitations of the potential-flow theory, of analyzing the 
flow in curved obstacles and generating a family of practi- 
cally suitable bucket shapes. Furthermore, this method has 
the inherent capability of being extended to the generation 
of even a larger family of concave-obstacle shapes through 
the use of distributed sources and sinks in the hodograph- 
plane . 

The method devised by Levi-Civita and successfully 
extended and employed herein requires the judicious selection 
of the variation of two parameters and the use of a suitable 
series (such as the Fourier series) to enable one to perform 
the integrals leading to the coordinates of the obstacle as 
well as of the free surfaces. A brief familiarity with this 
method allows one not only to generate a family of suitable 
concave obstacles but also to round the corners of a bucket 
composed of straight segments. This procedure has been 
illustrated by rounding the corners of a U-shaped obstacle 
previously analyzed by Sarpkaya [1]. 
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As to the Riemann-Hilbert mehtod, it has been demonstrated 
that it is in essence a mathematically more elegant and 
sophisticated version of Levi-Civita' s method in prescribing 
certain functions. It is for this reason that no particular 
solutions have been presented and only a linear variation of 
the parametric function was discussed. 

These methods have provided indirect solutions for the 
shape of the curved solid boundary in terms of a given jet- 
departure and jet-deflection angle. The direct ' solution of 
the corresponding problems requires the use of either a 
sequence of iterations or a numerical scheme for solving the 
Laplace-f ield problem. The latter method has been reserved 
for the analysis of the axisymmetric cases because of their 
far reaching and relative practical significance. 

The axisymmetric cases, namely, the deflection of a jet 
issuing from a nozzle from hemispherical thrust reversers 
have been solved directly through the use of the finite 
element method. In doing so, several new procedures have 
been introduced. These include the use of a self-adjusting 
and fixed grid and a convergent iteration scheme. The jet- 
deflection angle, the slope of the inner and outer surfaces 
of the deflected jet as well as the coordinates of the free 
streamlines were determined for a 'series of representative 
nozzle-hemisphere combinations. It has been demonstrated 
that the method and the iteration scheme are sufficiently 
general to analyze the more promising types of thrust-reversers , 
evaluate their performance characteristics, and to delineate 
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the range of importance of the geometrical variables involved. 
Finally, a few experiments have been carried out with air 
and water to assess the degree of significance of the real- 
fluid effects on jet deflection. These experiments as we 11 
as those conducted by others showed that in cases where the 
Coanda effect is negligible (water jet in air), the analytically- 
predicted and experimentally-measured jet-deflection angles 
are in good agreement and that the real-fluid effects are 
negligible. This is evidence of the fact that the deflection 
of a free jet by a solid boundary is well suited to potential 
flow analysis particularly for the cases where the inertia 
forces dominate the establishment of the flow pattern. In 
cases where the Coanda effect is not negligible (air jet in 
air) , the real-fluid effects are more pronounced and can cause 
the deflection angle to be as much as 20$ lower than those 
predicted from the potential-flow theory. Nevertheless, the 
potential-flow analysis of the two-dimensional as well as 
axisymmetric jet-deflection problems enables one to determine 
the optimum values of the geometrical variables involved and 
thereby eliminates the use of costly trial and error techniques. 
With such information available, the real-fluid effects could 
be evaluated with a few carefully chosen experiments and the 
refinements of design could be based upon a secure knowledge 
of the known principles of hydrodynamics. 
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c** 

C** APPENDIX A ** 

C ^ ^ 

C** COORDINATES OF THE BUCKET FOR K=.C0U86 CORRESPONDING ** 
C*»* TO A DEFLECTION ANGLE BETA=50.C CF.G. ** 



Q ^oi! :$ 4i £ 3 > * * r* Tjc 7 5 js « « A '■$ *■ *,'• * it £ Js w ti sf: ?.< i;; j£ :£ :>* ^ w -.$ :fe ajt jfc j; rf: * # j)t *> >? sj >jt >,t .* s< # 4 $ >> * 



X= 0 . 0 5C5 
X^C. 1002 
X = 0 .1531 
X=0 .2025 
X=0 .2577 
X=0 .3033 
X = 0 .3503 
X=0 .4008 
X=0 . 45 1 1 
X = 0 .5020 
X = 0 .5505 
X=0 .6026 
X = 0 . 6 5 3 9 
X = 0 .7061 
X = 0 .7513 
X = 0 .8048 
X=0 .3525 
X = C .9095 
X = 0 .9545 
X= l .0073 
X= 1 .0541 
X= 1 .1092 
X= 1.1523 
X=1 .2024 
X=1 .2622 
X= 1.3362 
X= 1 .3809 
X= 1 .4332 
X=L .4958 
X= 1 .5741 
X= 1 .6783 



Y=0. 7916 
Y=1 . 1097 
Y=1 .3654 
Y= l . 5637 
Y= 1 . 7553 

Y = l .8970 
Y=2 . 0305 
Y= 2 . 1626 

Y =2 . 2845 
Y=2 .3996 
Y=2 .5025 
Y=2. 6070 
Y =2. 7041 
Y=2 .7979 

Y = 2 . 87 53 

Y = 2 . 9629 
Y=3 . 0377 

Y = 3 . 1229 

Y = 3 . 1875 
Y=3. 2606 
Y=3. 3226 

Y = 3 . 39 3 1 
Y=3 . 4463 

Y =3 .5059 

Y = 3 . 5742 
Y=3 . 6546 
Y =3 .7011 
Y=2. 7533 
Y=3. 8128 
Y=3 .8823 
Y= 3.9634 
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C* *>« 
C* APPENCIX e *# 
c* ** 
C* MODIFIED HODOGRAPH METHOD ** 
E v 

C* COMPUTATION OF THE COORDINATES OF TFE BUCKET ** 
C* FOR A GIVEN DEFLECTION ANGLE 3ETA ** 
C* *v 



Q*# *****>j^**:*K*:fc***3:£y*>*jf: *#**££>,<*#**:************ * v *** *5* ** ’4 ** 

C 

IMPLICIT REAL*S( A-F t O-Z ) 

P 1 = 2. 1 A 1 5 92 60 0 
0 F I = 1. .O/PI 
BETA=0.0 
BETA=20 .0 
C 

C EVALUATE THE BUCKET COORDINATES EACF 10 DEGREES OF BETA 
C 

CC 1 1=1,7 
BET A=BETA+ 10.0 
BE=BETA*PI/180.0 
C 

C EVALUATE THE PARAMETER K (XK) FROM EC. 117) 

C 

XK = DEXP <-PI*DS IN (BE )/ ( 1 .O-DCOS (BE))) 

W R I T E ( 6 , 1 0 ) 

W R I T E { 6 , 11) BETA , XK 
• 1C FORMAT (5X»' 

11 FORMAT! 5X f 'DEFLEC. ANGL. BETA= ' , F5 .2 , « P AR AM .K= • , F 12 . 10 
0 ,/ ) 

c 

C DEFINE XO FROM EQ.<15) AND VARY T FROM -L TC C 
C 

XC = CPI-;' D L CG ( (I .0+XK)**2/ (4 .0*XK ) ) 

T = 1.00C01D0 
DT = 0 .00 0 IDO 
Y = C . 0 

DC 2 J = 1 1 1COOO 
T = T - D T 

Y=2.* (1 ./( T + XK) + 1 ./ (T-H ./XK)-2./ (T + l . ) )#DT/DLOG(T )+Y 
C 

C CALCULATE AND WRITE X AND Y COORDINATES 
C 

X=GPI* DLOG ( ( T + XK) *(T+ L.O/XK)/ ( T+ l . 0 ) ** 2 ) -X C 
WRITE (6, 12) X,Y 

12 FORMAT ( 5X, ' X=' ,F8. 3 ,3X , ' Y = ' ,F8.3 ) 

2 CONTINUE 

1 CONTINUE 
STOP 
- END 
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APPENDIX C 



EVALUATION OP THE COEFFICIENTS OF THE 
FOURIER SERIES REPRESENTING 6(a) 



Let e(a) be represented by the serie: 



6(a) = E a cos na 
n=0 n 



n = 0,1, 2, 3, . . . ( A-l ) 



According to Fig. 6, e(a) is defined by 



6(a) 

6(a) 

6(a) 

6(a) 

6(a) 

6(a) 



-IT 

IT 

2 



a 



°2 °1 



- 2 - 



°2 a l 



JL 

2 

H 

2 



JL 

2 



= 0 



rr -a , t 



°2 a l a 2~°l-' 



0 < a < a. 



< a < a 2 



°2 < ° < 2 



i < a < 0-a 2 



TT-02 < b < TT“0^ 



ir-a^ < a < -rr 



(A-2 ) 



Multiplying Eq. (A-l) by cos na and integrating over a period 
one has , 



TT 



IT 



/ 6(a) cos na da = a^ f cos na da 

“IT “IT 



(A-3) 



Now making use of the symmetry, one obtains 



2 u 

a = - / 0(a) cos na da 

n tt 0 



(A-4) 



97 



The above equation can be written as the sum of six integrals 
corresponding to each range defined by Eq. (A-2). Thus, one 
has , 



2 

a = — 
n ir 



f (-tt)cos na da + / j 

0 °i 



° 2~°1 



- 2 - 



°2 °1- 



cos na da 



2i 

2 



t r-a. 



+ / (- ^-)cos na da + f (- ^-)cos na da 

IT 

°2 2 



TT-a^ 

+ f 

tt— a_ 



TT 

2 



■rr-a^ i 



a 2~°l 



°2 °1 



TT 

cos no da + / 0 cos nada 

TT— O., 



} 



(A-5) 



After the integration one obtains. 



a = — ~ — [ (cos na 0 - cos na n ) - (-l) n • (cos na 0 

n d / \ d. j. d 

n 



cos na-^) ] 



(A-6) 



for n = 1,2,3,.. 



or 



n 2 (a 2 ~a 1 ) 



(cos na 2 - cos na^) 



( A— 7 ) 



for n = 1,3,5, 7 
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For n = 0 3 one has 



a 



o 



TT 



2 



( A— 8 ) 



as expected from the average value of the function 6(c). 
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C* APPENDIX D ** 

Qt)c 

C* LEVI-CIVITA' S METHOD ** 

Q «J- Uj -J . 

C* COM f-UTAT I ON CF THE COORD I AAT ES OF TPE BUCKET AND ** 

C* FREE STREAMLINES, AND THE DEFLECTION ANGLE BETA ** 

C* FOR GIVEN VALUES OF S I GMA I , S I GMA2 AND K ** 

C* ** 



C * :****#*****#£#* 5* #**#:**** * s)t j|: A ,< + r: 5 ;: •& ts r|c ■.% ,\- j{t & 

C 

C READ THE VALUES OF SIGNAL (SI), SIGPA2 ( S Z ) AND K (XK) 

C 

R E A D ( 5 , 1 ) S 1 , S2 , XK 
WRITE (6, 2) SL , S2 » XK 

1 FORMAT (2F10.0) 

2 FORMAT (5X, 1 SIGMA 1=' , F9.5,4X, ‘SIGMA 2=' , F9 . 5 , 5X , *K= ' , F8 . 
84) 

C 

C MAIN DO-LCCP TO INTEGRATE EQS . ( 40 ) » <4 1 ) ANC (42) 

C 

X = 0. 0 
V=C.O 

P I = 3 . 1415926 
FI2=PI/2.0 
S = P 1 2 

DS=S/3C0.0 
DC 3 1=1,300 
S=S-DS 
C 

C LOOP TO EVALUATE THE FUNCTION H(S),EQ.<38) 

C 

H = 0 .0 

DC 4 J = 1 , 200 1 t 2 
T = J 

H = H+ ( 2 . 0/ ( S2- S 1 ) ) * ( CO S ( T* S 2 ) -C C S ( T * S 1 ) ) * S I N ( 7* S ) / T** 2 
4 CONTINUE 
C 

C CEFINE FUNCTIONS OF THE INTEGRAND 
C 

F 1 = 2.0*XK*XK*SIN( S ) * ( 1 . 0+ S l N ( S ) ) / ( 1 . 0-XK*XMCCS ( S )**2 ) 
S M A = P 1 2 *( S -S 1 ) / ( S 2- S 1 ) 

C A = COS ( SPA ) 

S A = S I N ( SMA) 

C 

C SELECT VALUES OF THETA ACCORDING TO EQ.(39) 

C 

IFIS.GT.S2) C A= 0 . 0 
I F ( S . G T . S2 ) SA = 1 . 0 
IF(S.LT.SI) CA= 1 . 0 
IF(S.LT.Sl) S A=0 . 0 
C 

C INTEGRATE BY SIMPSGN’S RULE 

C 

X=X+EXP(-H)*F1*CA*DS/PI 
Y = Y + EXP (-H)-F l*SA*DS/P I 
W R I T E ( 6 , 6 ) S , X , Y 

fc FORMAT ( 5X , • SIGMA =‘ ,F8.4,5X, , X=',F7.3,3X, , Y= , ,F7.3) 

2 CONTINUE 
C 

C EVALUATION OF THE ANGLE BETA 
C 

C=( 1.0-SGRT( 1 .0-XK*XK) ) 

BETA=-ATAN(2.0*C/( 1.0-C*C) ) 

DC 7 L = 1 , 301 , 2 
T= L' 

BETA= ( 2 ./( S2-SI ) ) * ( CO S ( T* S2 ) -C 0 S ( T*S L ) ) *C**L/T**2 + BETA 

7 CONTINUE 

BETA=18C. 0*( 1 .0+ BETA/ PI ) 

WRITE(6,8)BETA 

8 FCRMAT(5X» ' BET A= ' » F 10 .3 ) 



100 



non odd 



C EVALUATION OF THE C0ORDENATES OF THE FREE STREAMLINES 
C 

I NC = 0 
R=-0. OOCL 
CR=C/I000 .0 

12 CONTINUE 
X = C .0 

Y = C . 0 

DO II 1=1,1000 

p = R -f OR 

Fl= (XK-XKV (R”-R-l . 0/R**2 )/ ( 2.0*R ))/(( XK**2/^.0)^(R+ 1.0/ 
8R ) #*2-1 .0) 

EVALUATE W(R) FROM EG. (29) 

h=-*ATAN (2.0*R/( L.0-R*R) ) 

DC 10 J=lt 1001,2 
T = J 

W=W+(2.0/I S2-S1) )’MC0S(T*S2J-C0S(T*Sl) ) *R*» J/ T*«2 
1C CONTINUE 

X=C0S( W }*F1*DR/P I+X 

Y = S 1 N ( W )*F1*DR/P I + Y 
WRITE {o,I3J R,X,Y 

13 FORMAT ( 8X . 'R=* ,F 10.6, 5X, 'X= ' ,F9.3,2X, ' Y= • , F9.3 ) 

11 CONTINUE 

GO TO 12 TO EVALUATE THE UPPER FREE STREAMLINE 

I ND= I ND+ 1 
R= 1 .000 1 

0R= (C-l .01/1000.0 
I F ( IND. EQ. 1) GO TO 12 
STOP 
END 
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APPENDIX E 



EVALUATION OP THE INTEGRAL FOR w(t) 



w(t ) 




k 

/ 

1 



-(a + bn) 
(n-t) \/^n(n-k) 



dn 



(E-l ) 



Noting that 



a , bn a + bt , , 

— + ~r - z — + b 

n-t n-t n-t 



Eq . (E-l) can be written as. 



co(t ) = 



- V^t (t-k) 



IT 



k 



(a + bt ) / 



dn 



+ b / 



d 



l (n-t) |/- n(n-k) l / -n(n-k) 



(E-2) 



The second integral in Eq. (E-2) is 



k 

f 

1 



dn 

]j n(k-n) 



cos 



-1 2 - k 
k 



(E-3) 



The first integral in Eq. (E-2) depends on the particular range 
of t. For 0 < t < k , one has. 
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1 



In 



dn 



| (t+1) - t + ^k-1 \[t{ k-t) 



l (n-t) j/ n (k-n ) ^ t(k-t) 



| (t-i) 



(e-*d 



For t < 0 and t > k one obtains. 



dn 



l (n-t) \j n(k-n) ^Tt(k-t) 



cos 



2t 
-1 k 



- t + l 



t - 1 



(E-5) 



Inserting Eqs. (E-3) > (E-*l) and (E-5) in Eq. (E-2) one finally 
obtains , 



k 

/ \ a + bt , 2 

w(t ) = In — 



(t+l) - t +\/kTVt(k-t) b j — ■ 2 - k 

y t(k-t) cos — j— 



- 1 (t - 1) 



for 0 < t < k, and 



2t 



i(a + bt) -lk t - 1 £ ^ ~ -1 

co(t) = — — - — - cos — ut(t-k) cos 



-1 2-k 



IT 



t-1 



for t < 0 and t > k . 
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APPENDIX F 



EVALUATION OF THE DEFLECTION ANGLE 3 FROM EQ. (52) 



Letting t - vco in Eq . (52) one has 



lim w(t) = 



i / (at bn) . dn 



11 1 \j ri ( k-ri ) 



(F-l) 



and noting that 



3 = Tr-lim iw(t) 
t-t-co 



( F-2 ) 



Eq. (F-l) can be rewritten after some rearrangement as. 



3 = it + 



/ . bk N 

( a + ~2*) 

IT 



k 



/ 

1 



dn 



(k-n ) 



7 ,/n(k-n) 



k 



(F-3) 



Evaluating the integrals and applying the limits of integration, 
one has, 



3 = it + — 

IT 



, , bk N -1 2 - k , , f, 7 

(a+-p-) cos — j- — + b^k-1 



(F-'O 



Finally, replacing the corresponding values of a and b one 
obtains , 



3 



IT + 



1 

2 (k - 1) 




2 ) 



cos 



- 1 2 - k 
k 




10H 



APPENDIX G 



LIST OP COEFFICIENTS K e . AND B e . 

ij J 



s ii ■ (a i a j + b iV /60Ae 



*h- 


' 3S n 


(3^1 + r 2 + 


r 3> 














K 12 ' 


■ K !l 


S 12 ( ‘ 2r i 


+ 2r 2 


+ r 


3 } 










k 13 


‘ K 31 


= - s i 3 ( 2r i 


+ r 2 


+ 2r 


3 } 










k i4 = 


■ 


s 11 (3r 1 


2r 2 - 


r 3 ) 


+ 


S 12 


(l^r 1 


+ ~ >r 2 


+ 3r 3 ) 


K 15 ' 


■ K 51 


' S 12 (3r l - 


r 2 


2r 3 ) 


+ 


S 13 


(3r x - 


■ 2r 2 ■ 


‘ r 3> 


K l6 = 


= K 61 


' S ll (3r l - 


r 2 


2r 3 ) 


+ 


S 13 


( l^r ]L 


+ 3r 2 


+ 3r 3> 



K 



22 



3S 22^ r l + 3r 2 + 



K 23 


= k 32 


= - s 23 ( r]L + 


2r 2 + 


2r 3 ) 










K 24 


= k 4 2 


■ S 12 (3r l + 


li)r 2 + 


3r 3 ) 


+ 


s 22 ( 


2r i 


+ 3r 2 - r ) 


LT\ 
CD OJ 


K 52 


" S 22 ( - r l + 


3r 2 - 


2r 3 ) 


+ 


S 23 (3r 


i + 


l^tr 2 + 3r 3 ) 


K 26 


= K 62 


= S 12 (- r x + 


3r 2 - 


2r 3 ) 


+ 


S 2 3 ( _ 


2r i 


+ 3r 2 - r 3 ) 



K 



33 



3S^2( r 1 + r 2 + 3r^) 



K~ i, = K*L = S 10 (- r, - 2r n + 3r 0 ) + S 00 (- 2r, - r 0 + 3r,) 



‘3^ 



K 



l 43 13 A 1 “2 ’ Ji 3' ' °23 v “1 "2 ' Ji 3’ 

K® = S 00 (3r n + 3r 0 + l^rj + S„(- r, - 2r 0 + 3rJ 



35 53 23 



33 



K e 36 = Kg 3 = S 13 (3r 1 + 3r 0 + l4r 0 ) + S 00 (- 2r n - r 0 + 3r„) 



33 
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k 44 = 8 ^ S ll^ r l + 3r 2 + r 3^ + S 12^ 2r l + 2r 2 + r 3^ 

+ S 22^ 3r ] + r 2 + r 3^ 

k ^ 5 - K 54 - 8s i 3 < r i + 3r 2 + r 3 > - *1(3^ + S 22 r 2 + S.,^) 

K l|6 = K 64 = 88 23^ 8r i + r 2 + r 3^ " l| ^ s ll r l + S 12 r 2 + S 13*V 

K®g = SESjjfrj^ + r 2 + 3r ) + S 23 ^ r i + 2r 2 + 2r 3 > 

+ S 33 (r 1 + 3r 2 + r 3 ) 

K^6 ~ “ 8S-^2 ( r -L ^ r 2 * 3r 3^ ~ ^ ^ S i3 r i "*" ^23^2 ^33^3^ 

K 66 = 8 ^- S ll^ r l + r 2 + 3l V + S 13^ 2r l + r 2 + 2r V 

+ + r> 2 + r^) 



In the matrix B® terms, 2,^ = length of side i of a triangular 



element. Thus 



B 



e 

1 




+ 





1 

z 



B 



e 

2 






+ 





1 

z 



B 



e 

3 






1 

z 
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B 4 ^ r l + r 2^9n^ ^3 3 



B e = (r + r ) (■^•) — 

5 ^2 3 U 3n ; 1 3 



B e = (r + r )(— ) a2 £ — 

B 6 ^ 3 l n 3n ; 2 3 
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£***$*£ £:£# i<% ## v* 



C** * * 

C** APPENDIX H ** 

C*# A V 

C** CGORD I NAT FS OF FREE- STREAMLINES ** 

>>* 



Q^^********************'************** ************** ********* 



CASE NUMBER ONE 



BE TA = 20. 3 


DEG 


V0=0. 817 




X=- 0 . 400 


R=1 .000 


X=-0. 370 


R=1 .003 


X=-0 . 325 


R= 1 .009 


X=-0. 301 


R=1 .012 


X=-0 .263 


R= 1.020 


X = -0 .230 


R=1 .030 


X=-0. 203 


R= L . C39 


X= — 0 . 181 


R=1 .042 


X = -0. 130 


R=1 .060 


X=-0.C99 


R= l .080 


X=-0 .066 


R=1 .100 


X=-0. C38 


R= 1 » 1 1 8 


X=-0 .016 


R=L. 137 


X= 0.001 


R=1 . 158 


X= 0.018 


R=l. 176 


X= 0.034 


R= 1 . 192 


X = 0.C53 


R=1 .224 


X= 0.063 


R= L .258 


X= 0.071 


R = 1 .290 


X= 0.C72 


R= 1 . 323 


X= 0.072 


R = 1 .340 


X= 0.C69 


R = 1 .357 


X= 0.C65 


R= 1 . 37 5 


X= 0.0 59 


R= 1 .392 


X= 0.056 


R= 1 .398 


X= 0.053 


R= 1.405 


X= 0.050 


R= 1 .412 


X= 0.C47 


R=1 .419 


X= 0.044 


R= 1.425 


X= 0.040 


R=1 .431 


X= 0.036 


R= 1.438 


X= 0.032 


R= 1 .444 


X= 0.028 


R=1 .450 


X= 0.023 


R = 1 . 4 5 6 


X= 0.018 


R= 1 .462 


X= 0.C13 


R= 1 .463 


X= 0.007 


R= 1.474 


X= 0.001 


R = 1 .48 1 


X=— 0 «• 006 


R=l .480 


X= — 0 . 0 l s> 


R= 1 .495 


X=-0. C24 


R= 1 .505 


X=— 0 . 036 


R= 1.514 


X=-0 .048 


R= 1.524 


X=-0. 062 


R= 1 . 533 


X=-0 .076 


R= 1 . 544 


X=-0 .095 


R= 1. .55 5 


X=-0 . 1 13 


R= l . 565 


X=-0 . 131 


R= 1.576 


X=-0. 183 


R=1 .607 


X=-0 . 252 


R= 1.634 


X=-0 .328 


R=L .669 


X=- 0.385 


R=1 .694 


X=-0 .457 


R= 1.724 


X=-0 .530 


R=1 .753 


X=- 0 . 6C0 


R= i . 7 8 1 


X=-0 .658 


R= 1 .805 
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X= — 0 .014 R= 1.805 

X =- 0 . C 3 0 R= L .80 7 

X=— 0 . C4 5 R = L . 3 1. 3 

X=-0 .061 R = 1 .315 
X=-0 .111 R= 1.834 

X=— 0.173 R= 1.855 

X=-0 .252 R= L . 884 
X=-0 « 310 R= 1 . 906 
X=-0 . 373 R= L . 933 
X = ~0 . 457 R=1.959 

X=- 0.52 8 R= 1.934 

X= — 0 .58 7 R= 2 . 00 5 



CASE NUMBER TWO 
BETA = 3 5. 0 DEG 



V0=0. 885 

X=-C.800 
X=-0 « 7 74 
X=— 0 .743 
X=-0. 723 
X=-0 .697 
X=-0 .663 
X=-0 .629 
X=-0 . 595 
X=-0. 56L 
X=- 0 .525 
X=-0 .489 
X=-0. 453 
X= - 0 . 4 1 8 
X = -0 .394 
X=-0. 270 
X= — 0 .346 
X=- 0 . 324 
X=-0.30i 
X=-0 .281 
X=— 0 . 260 
X=-0 .24 1 
X=-0 . 224 
X=-0. 208 
X=-0 . 194 
X=-0 . 181 
X=-0 . 166 
X = -0 .155 
X=-0. 146 
X=-0 . L42 
X = -0 . 140 
X=— 0 . 140 
X=-0 . 143 
X=-0. 149 
X=- 0 . 154 
X = -0 .161 
X=— 0 . L68 
X=— 0 .177 
X = -0 . 184 
X=-0. 191 
X= — 0 . 200 
X=-0. 209 
X=-0. 243 
X = -0 .304 
X=— 0 .367 
X=- 0.430 
X=-0 .494 
X=- 0.562 

X=-0 .CCO 
X=-C .068 
X=-0 . 083 
X = -0. 137 
X=— 0 . 150 



R=1 .000 
R=1 . 001 
R= 1 .003 
R= l . 006 
R=1 .009 
R=1 .013 
R= 1 . C 1 9 
R= 1.025 
R=i .030 
R= 1.038 
R= L .047 
R= 1 . 057 
k= 1 . 069 
R=1 .080 
R=l . C89 
R= 1.100 
R = 1 . 11 1 
R= 1.12 5 
R= 1 .139 
P.= 1.153 
R= 1.168 
R=L . 182 
R=1 . 199 
R= l .216 
R=1 .235 
R- 1.258 
R= 1 .282 
R=1 .308 
R=1 .335 
R=l .358 
R= 1.382 
R= 1.404 
R=1 .425 
R= 1.443 
R= 1 .461 
R=1 .47 8 
R= 1.495 
R=1 .506 
R = l . 5 1 7 
R=1 .529 
R=l .539 
R= 1.579 
R= 1 .635 
R=1 .689 
R=1 .733 
R=1 .79 1 
R= l . 838 

R= 1 .800 
R=1 . 8C0 
R= 1.803 
R=1 .809 
R= 1 . 855 
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X=-0 . 234 R= 1 .88 2 
X=- 0.208 R= l .947 
X = -0 .26 7 R = 1 . 9 7 6 

X=-0 . 425 R = 2 .03 7 



CASE NUMEER THREE 



BETA= 50 . 5 
V0=0. 810 

X=- 0.800 
X=-0 . 765 
X=-0. 729 
X=-0 .694 
X=— 0.660 
X=- 0.632 
X=-0 .603 
X=-0 .575 
X=-0 . 548 
X=-0 .525 
X=-0. 5C3 
X=- 0 .482 
X=-0 .461 
X=-0. 451 
X=-C .440 
X=-0.430 
X=— 0.420 
X = -0 .4 12 
X=-0.4C4 
X=-0 .296 
X = -0 .380 
X=- 0.384 
X=-0 .379 
X=-0 .375 
X=— 0.37 l 
X = -0 .366 
X=-0. 361 
X=— 0 .257 
X=-0 .353 
X=— 0.349 
X=-0 .346 
X=-0 .343 
X=— 0 . 24 L 
X=-0 .339 
X = -0 .337 
X=~0 .235 
X=-0 . j34 
X=— 0 .333 
X=-0 .232 
X=-0 .331 
X=-0 .221 
X=-0 .327 
X=-0.367 
X=- 0 . 399 
X=-0 .495 
X=-0. 561 
X=-0 .627 

X=-0. 000 
X=-C .C90 
X=-0 . 101 
X=- 0 . 158 
X=— 0 .175 
X = -0 . 25 5 
X=-0 . 2 22 
X= — 0 .397 
X=-0 .469 



DEG 



R=1 .000 
R= 1 .007 
R=i .015 
R= 1.025 
R = 1 .036 
R= 1.045 
R= 1.055 
R= 1 . 0t>6 
R= 1.07b 
R= 1 .090 
R = 1 .102 
R= l . 1 15 
R=1 .130 
R= l . 138 
R= 1 . 147 
R=L . 155 
R= 1 . 165 
R= 1.17 3 
R=i . 182 
R=l. 191 
R=1 .200 
R=1 .206 
R= 1 .212 
R= 1.219 
R- 1.225 
R= 1.234 
R=1 .242 
R- 1.257 
R= 1 .260 
R=1 .268 
R= 1.276 
R=l .204 
R=i .293 
R= 1.300 
R= L .308 
R= 1.316 
R = 1 .325 
R= 1 .331 
R- 1.337 
R=1 .344 
R=i .350 
R= 1.443 
R= 1 . 5 1 3 
R= 1 .61 1 
R=1 .742 
R=1 .830 
R=1 .910 

R=1 .600 
R=1 .580 
R = 1 .654 
R= 1 . 67 8 
R= 1.734 
R=1 .809 
R= 1 . 697 
R= 1.977 
R = 2 . 059 
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^ **£# #54#* >* # •£ * sjsjic A # # 2s# * # 



C ^ ^c^c 

C* APPENDIX I ** 

C* >*f 

C* FINITE ELEMENT METHOD 

c ^ ^ # 

CM COMPUTATION OF THE ANGLE CF DEFLECTION AND TEE M* 

CM COORDINATES OF THE FREE STREAMLINES FOR AN ** 

CM HEMISPHERICAL BUCKET ** 

CM s* 



C M M M M M '.<■ 4 >? M M •!( v A is •* ;|s ^ £ v & ^ # £ # # j}e is £ *; sfc * # Hi £ % * # * R # s-s :£ >, L ^ $ :)< ^ # :)c # ^ >* *: # 

c 

IMPLICIT REAL*8< A-H,0-Z) 

C 

DIMENSION IN0(288,6),JET(9,3),JABA(12,6),X(683),R(682) 
8,RR(65),XX(65) f ELE(6,6),Z(683,60) ,B(683) , ANGLE (45) 

8 , V ( 80 I 
C 

C VALUES FOR THE NUMBER OF ELEMENTS ( N'UMEL J » 8 ANDh IDTH ( N BAN) 
C NUMBER OF NODAL POINTS (NNODE ) , NODE AT LIP OF E LCKET (IM2) 
NNCDE = 6 83 
N E AN=60 
N2=597 
NLMEL=2S6 
PI=3. 14 L5926D0 
DK=1. ODC 
C 

C ASSIGN STARTING VALUES FOR DEFLECTION ANGLE (BETA) AND 
C VELOCITY AT THE ENTRANCE OF THE NOZZLE (VO) 

BETA=0. SCDC 
VC=0 . 85 CO 
F ACU=0 . 001 DO 
FACL=0 » C2D0 



C 

C READ ECU I VALENCE CF NODAL POINTS 

RE AD ( 5 , 1 ) ( I ,( INO( I , J > , J = 1 , 6 ) , N = 1 , NUME L J 

1 FORMAT (3(713 ) ) 

C 

C READ COORDINATES OF INTERIOR POINTS 
RE AC (5, 2) ( I,X( I),R{ I) ,N=1,90) 

2 FORMAT! 5( 13, F6.0-F5.0) ) 

C 

C READ X VALUES OF THE NODES OVER THE BUCKET 
READ( 5,3) ( I , X ( I ) , N = 1 ,18) 

3 FORMAT ( 9( 1 3, F5 « 0 ) ) 

C 

C READ EQUIVALENCE CF NODES IN THE JET SECTION 
RE AD (5, 4)- ( ( JET( I , J ) , J=l,3) , 1 = 1 ,9) 

4 FORMAT (8(313) ) 

C 

C READ ECUIVALENCE OF NODES IN THEFLOATING GRID 
R E AC ( 5 ? 5 ) ( <JABA(1,J) , J = 1 , 6 ) , 1 = 1 , 1 2 ) 

5 FORMAT (4(613)) 

C 

C. READ ThE COORDINATES CF THE ASSUMED FREE SURFACE 
R E AD ( 5 , 6 ) (I ,X( I ) ,R( I ) ,N=1, 10) 

6 FORMAT (5(I3,F6.0,F5.0) ) 

C 

C LOCATE THE NODES OVER THE BUCKET 
R AC I 0= L .6000 
DC 40 1=63,78 
J= INC ( I ,3) 

4 0 R( J)=USCRT(RADI0**2-X( J)”-*2 ) 



C 

C 



LOCATE THE UPPER FREE SURFACE REALTIVE TO POINT N2 
PX=0.077D0 
D I SP=0 . 0 
DC 10 1=55,59 
XC=6I— I 
K= I NO ( I , 1 ) 

X (K)=X (N2 l-CXMXD 

IC R(K)=R( N2) -DTAN(BETA)* ( X ( K ) -X ( N2 ) + C I S P ) 



111 



c 

C LCCATE THE PCIMS CLOSE TO THE LIP 
ZP=R(G47)-R(N2) 

X (617 )=X(647) ^0. 2500 

X( 627) = X( 6VD*0.5D0 

X<637 ) = X( 647DG.75D0 

F(617)=R(N2)+ZR*0.25C0 

R ( 627 ) = R ( N 2 ) + ZR* 0 . 50D0 

R (637 ) = R ( N2 ) + ZR*0 .7 500 

R E A U ( 5 , 6 ) { I ? X ( I } , R ( I ) , N= 1 , 1 6 ) 

NAL = 9 
C 

C STARTING POINT FOR ITERATIONS ON (BETA) AND (VC) 

10CC CONTINUE 
INC EX = 0 
C 

C LOCATE THE REGION C OF THE LONER FREE SURFACE 
T R E S = 0 . 05D0 
VETA=BE TA + TRE S 
DC 20 1=1, MAL 
1 1 = JET I 1,1 ) 

I 2 = JET( 1,2) 

I 2 = JET ( 1,2) 

I F ( I ,6c .8 ) GO TO 21 

R ( I2)=DSGRT(R(I 1 ) **2-VC*DC0 S ( VE T A ) ) 

X ( 1 2 ) = X ( I 1 ) - DT AN ( V E T A ) # ( k ( I 1 ) ~ R ( 12) ) 

21 X(I3)=(X(J1)+X< 12) ) *0.5 CO 

2 C R(I2)=(R(I1)+R( 12) DC.5D0 

X (5S1 ) = 0 . 06 DO 
R (551 ) =1.4800 
X ( 606 ) =-0 . 16D0 
P (606 ) = 1 .4800 

I F ( MAL . EQ. 9 ) X(618) =-0.3300 
I F { MAL . EG. 9 ) R (6 18 )= 1.44D0 
C 

C STARTING POINT TO ITERATE THE FREE SURFACES 
100 CONTINUE 

I N C E X = I NDEX+1 
C 

C LCCATE THE INTERIOR POINTS IN THE F L CAT I NG GR1C 
OC 30 1=2,12 
I 1 = JABA ( I , 1) 

I 2 = JABA ( I ,2) 

I 2 = JABA ( I ,3) 

I 4 = JABA ( I ,4) 

I 5 = JABA ( I ,5) 

1 6 = JA B A (1.6) 

X( I2)=X( I I ) + 0.l5DO#(X( I4)-X( I 1) ) 

X ( 1 3 ) = X ( I 1 ) +0 . 5 DO* ( X ( 14 )-X (ID) 

R( 1 2 ) = R ( I 1) + 0.15D0-(R( I 4 ) - R ( I 1) ) 

R ( 13) =R ( I 1 M0.50D0*(R( I4)-R( ID ) 

X ( I 5) =0 .500* ( X( I 6 ) +X ( 12) ) 

3 C R ( 1 5 ) = C . 5 D C* ( R ( I 6 ) + R ( 12) ) 

C 

C LCCATE INTERNE Cl ATE PCINTS IN ASSUMED LOWER FREE SURFACE 
DC 50 1 = 1, 1.2 
K= J ABA (1,1) 

X X ( I ) = X ( K ) 

5C R R ( I ) =R ( K ) 

CALL DRVF' IT(2, 12, RR ,XX , 4 , 4 5 , RR , XX, ANGLE ) 

DC 60 1=9,48 
K = INO ( 1,3) 

X (K)=XX (1-3) 

6 C R ( K ) =RR ( I - 3 ) 

IF (KAL .NE..5) GO TO 28 
X (618) =-0.3300 
R(618) = 1.4 4D0 
X ( 628 ) =-0 . 36D0 
R (628) =1 . 52D0 
X(638)=-0. 4100 
R ( 63 8 ) = 1 . 6 ODO 
28 CONTINUE 
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RAL' = 5 7.295800 

ALFA3=DATAN( ( R( 6 72 ) -R ( 647 ) ) / ( X ( 64 7 )-X ( 672 ) ) ) * R AD 
W R I T E ( 6 , 6 7 ) ALFA3 , ALFA3 
C CALCULATION OF THE DEFLECTION ANGLE ALFA 

ALF A 1= DAT AN ( ( R ( 6 72 )-R ( 667 ) ) / ( X ( 66 7 ) -X ( 6 72 ) ) ) * R AD 
ALFA2=DAT AM ( R(6 77 } -R ( 6 57 ) ) / ( X ( 6 5 7 )-X ( 6 7 7 ) ) )*RAD 
GAMA1 = DATAM ( R I 6 6 8 ) - k < 64tt ) ) / I X I 64 8 ) --X 1 66 8 ))) ’> R AD 
GAFA2=DATAN( I R I 6 5 8 ) ~R I 638 i ) / 1 X I 63 £ )-X I 65 8 ) ) )*RAD 
l\RITE(6,6L ) A L F A L , G A M A 1 
WRITE(6,61) ALFA2 » GAMA2 

61 FORMAT ( 5X * * UPPER f SL * , F 1 0 . 2 , 5X , • LOWER FSLSF10.2) 

G AM A3 = DATAM I R I 66 6 ) -R I 658 ) ) / ( X { 6 58 ) -X ( 668 ) ) )*RAD 
G AM A4=DAT AM I R I 65 8 ) -R ( 648 ) ) / ( X l 64 8 )-X ( 6 5 8 ) ) ) *P AD 
WRITE (6, 67) GAMA3 , G AMA4 
67 FCRMAT(2F2G.2) 

C 

WFITE(6»500) VO, BETA 

5CC FORMAT ( / , 5X » ' V0 = ’ , F 1 C . 4 , 5X , 1 BE T A= • , FI 0 . 4 , / J 
C 

C ZERO CUT MATRICES (Z) AND (B) 

DC 42 L=1 iNNODE 
DC 43 M--1.NBAK 
43 Z I L , M ) = 0 . OCO 
4 2 BIL) = C . CDO 
C 

C EVALUATE THE ELEMENTS CF MATRIX (Z) 

DC 44 K = 1 , N U M E L 
1 1= INO ( K , ]. ) 

12=1 N OIK, 2) 

I 2 = 1 NO ( K , 3 ) 

A 1 =X I I3)-X( 12) 

A2 = X( 11 )-X( 13) 

A 2 = X( I 2 ) - X I I 1 ) 

E L = R I 12 )-R ( 13 ) 

B2 = R( I 3 ) - R I I 1 ) 

B 2 = R I I1KRII2) 

AREA=0.5*( A2*B1-A1»B2 ) 

SI L=( Al*Ai+Bl*Bl )/ I AREA*60.0) 

SI 2= I A1*A2+BL*B2) / <AREA*60. 0) 

S 13 = 1 A1*A3 + 81*B3 )/ (AREA*60.0) 

S22=( A 2*A 2 -f-b 2* B2 ) / (A REA ”60.0 
S23=(A2*A3+B2*63)/ (AREA*60. C) 

S33 = ( A3*A3+B3*B3 )/ (ARE A *60.0) 

ELEI 1, 1 )=3.0*Sli*( 3.0*R(I IJ+R(I2)+R(I3)) 

ELEI L,2) = ~SL2 y< (2.C*R( 1 1 J+2.0*R( I2)+R( 13) ) 

ELEI2, 1 )=ELE(1 ,2) 

ELEI L , 3 ) =•— S 13" ( 2 . C * R I I l)+R( I 2 ) + 2 . C*R I 1 3 ) ) 

ELEI3, I )=ELE( 1,3) 

ELE(L,4) = S11"(3. 0 >; R ( I 1 )— 2.^RI I2J—RCI3) )+S12^I14.0 s f‘R( 
8 I l )+3.0*R( 12 )+3.0*R( 13)) 

ELEI 4, L ) = E L E I 1 , 4 ) 

ELEI 1 , 5)=S12M3. -'R I I 1) - R< 12 >-2.*R( 13) ) +S 13 » 1 3 . 0*R I I 1 
8 )-2 . 0 R 1 12 >-RI 13 ) ) 

ELE(5,1)=ELE(1,5) 

ELEI l,6)=Sil*(3. 0*R( II )-R( 1 2 ) -2 .C*R( 13) ) + S13-M 14.0*R 
8(I1)+3.C*R(I2)+3.0<R(I3)) 

ELL(6,1)=ELE( 1,6) 

ELEI 2, 2 ) = 3.C*S22M R( I 1)+3.0*R( I2)+R( 13) ) 

ELEI2 ,3 )=-S23*< R( I 1 )+2 .0*R( I2)+2,0*R(I3)) 
ELE(3,2)=ELE(2,3) 

ELEI 2 * 4 ) = S12*(3.0*R( I 1)+14.0*R( I 2 ) +3 . 0*R I I 2 ) ) + S22* I- 
8 2 . 0 J : R I I 1I+3.0VRI I2J-RI 13) ) 

ELE(4,2)=ELE(2,4) 

ELEI 2, 5 ) = S22*(-R( I 1)+3.0*R( I2)-2.0*RI 13) ) + S 23* I 3 . *R I 
8 I i ) + 1 4 . C* R I 1 2 ) +3 . 0 * R I I 3 ) ) 

ELE(5,2)=FLE(2,5) 

ELE (2 , 6 )=S12*(-R( I 1 H3.0~<R( I2)-2.C*R( 13) )+S23*(-2 . C* 
8RI 1 1 ) +3 . C”R I 12J-RI 13) ) 

ELE(6»2)=ELE(2,6) 

E L E 1 3 , ) = 3 . X: S33' : I R I lll+RI I 2 ) +3 . 0~R ( 13) ) 

ELE(3,4)=S13*( -R ( I 1) -2.0*R( 12) +3 .0*R I 13 ) )+S23*I~2 .0* 
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on on oooo 



SR ( 11 )-R ( 12 )+3 .0*R ( 13) ) 

ELE( 4, 3 )=ELD <3 ,4 ) 

ELE (3, 5)=S23*<3 .0*R ( 1 L)+3.0*R< I 2) + 14 . 0*R ( I 2 ) 1+S33* (- 
8R (11) — 2 e 0 M I2)+3.0*R( 13) ) 

ELE( 5,3)=f:LF (3,5) 

ELE(3,6)=S13*(3.0*R( I i)+3.0*R( 12)+14.0*R( I2))+S33*(- 
82 . 0*R ( ) 1)-R( 12) +3 ,0*R( 13) ) 

ELE(6,3)=ELE(3,6) 

ELE(4,4)=8.0*( Sll*(R( I 1 )+3.0*R( I2)+R( 13) ) + S 1 2* ( 2 . 0*R 
8 ( I 11+2 . C*R( 12 ) +R( 13) ) +S22--M3 ,*R ( 1 i )+k( 1 2 ) 4p ( 13 ) ) ) 

ELE (4, 51 = 8.0*S13-M R( I 1 ) +3 . 0*R ( 1 2 ) + R ( 1 3 ) ) - 4 . * { S 12*R ( 1 
81 ) + S22'<R( I2)+S23*R( 13) ) 

ELE(5,4)=ELE(4,5) 

ELE (4, 6 )=B ,0*S23*( 3.0*R( 1 L ) +R ( 1 2 i +R ( 1 3 ) ) -4 . * ( S 11* R ( 1 
81) + S12-R( I2)+S13*R ( 13 ) ) 

ELE(6,4)=ELE:(4 r 6) 

E L E ( 5 , 5 ) = 8 . 0 'M 5 2 2” ( R ( I 11 +R ( I 2 ) + 3 . C*R ( 1 3 ) ) 4 5 23* ( R ( I L ) 
8+2.0*R( I2)+2.0*R( I3)+S33*(R( 1 1 ) +3 .*R< I 2 ) +R ( 13) ) ) 

ELE( 5, 61=8.0*5124 ( k( I 1)+R( I 2) +3. C*R( 13) ) -4 . C* ( S 13* R { 
8Ii+S23*R(I2)+S33*R(13) i 
ELE(6,5)=ELE(5>6) 

ELE(6, 61=8.0* (SI i*(Rl I 1 } +R { I 2 J + 2 . C*R ( 13) )4S13*(2.0*R 
8 ( 11J + R( 12) *2.*R( 13 > )+S33*(3.*R( 1 1 ) +R ( 1 2 ) +R ( 13 ) ) ) 

C 

C L C AD MATRIX (Z) 

DC 45 J = 1 , 6 
I I = I NO ( K , J ) 

DC 46 L = 1 , 6 
I I 1 = 1 NO (K, L ) 

I F l I I.GT.I II } GO TC 46 
NN* 1 1 1- 1 1+ 1 

IF(NW.GT.NPAN) NBAN=NW 
Z(II,NW)=Z(II,NW)+ELE(J»L) 

4 6 CCNTINUE 
45 CCNTINUE 
44 CONTINUE 

LCAD VECTOR ( B ) KITH BOUNDARY CONDITIONS 

AT THE INLET OF THE NOZZLE 
B ( 1 ) =— 0 . 0 

B(2)=-R(3)*R(3)*VO/3.C 
B(31=-R(3)*R(5)«V0/6c0 
B ( 4 ) =- ( Rl 3 ) +k ( 3 ) ) * ( R ( 5 ) -R ( 3 ) ) ,: VC/3.0 
B ( 5 )=tR ( 5 ) M k ( 7 ) — R ( 3 ) ) *V0/6 . G 
B(6)=-(R(7 HR(5) )*(k(7)-R(5) )*VG/2.0 
B(7)--R(7)' : (R(9)-R (5) )” V0/6.0 
B(8)=~ ( k(9)+R(7 ) ) v (R(9 )-R ( 7 ) ) "VO/3 .0 
6 ( 5 ) =-R (9 ) * ( R ( 9 ) -R i 7 i ) * VO/ 6 .0 



AT TEE EXIT OF THE JET (TAKE POTENT I AL = 100 . 0 ) 

DC 70 1=668,672 
B ( 1 ) = 10,0D12 
70 Z ( 1 , L)= 10.GDL L 

SOLUTICN OF THE POTENTIAL BY GAUSSIAN ELIMINATION 
M A L = 2 

CALL BANDEC( Z ,B, NNCCE, NBAN, 1,683,60,1) 

C 

C CALCULATIONS FOR THE VELOCITIES ALONG THE BOUNDARIES 
DC 110 1=1,79 
J 1= INO ( 1 , 1 ) 

J3 = I NO ( 1,3) 

S = D SQR T ( ( R ( J 1 ) -R ( J3 ) ) **2+ ( X ( J 1 ) -X ( J3 ) )**2) 

V( 1 ) = ( B( J3)-B( J 1) )/S 
I F ( I . GE . 55 ) V ( I ) =- V ( 1 ) 

IF(I.GF.55) J3 = J 1 

to F I T E ( 6 » 111) J3iX(J3),R(J3) , V ( I ) 

111 FCPPAT (5X» ‘NODE' , 15 , 5X » * X=', F 7 .3 , 3X , • R = • , F7.3,10X, 
8 V ELOC= ' , F 7 . 3 , / ) 

110 CCNTINUE 



C VELOC. AT LIP CF NOZZLE (VI) AND LIP OF BUCKET (V2) 
V1=(B( 194J-B! 199) )/ 0.0 A 75 DO 
V2=V( 63 ) 

WRITE<6,112) VI, V2 

112 FORMAT!//, ' VELOC. L I P 1 = ' , F 1 0 . 4 , 5X , ' VELCC. L I P 2 = ' , 
6F1C.4,//) 

C 

C CORRECTION FOR THE UPPER AND LOWER FREE SURFACES 
WRITE (6, 122) 

C 

C FCP THE LOWER FREE SURFACE 
DC 120 1=9,54 
AN=PI-dETA 

1FII.LE.48) AN=ANGLE( 1-4) 

K = I NO ( I » 3 ) 

DELV= ( V (I ) **3- 1.0) *FACL 
WRITE! 6, 121) K , U E L V 
XXII )=X(K 5 -DEL V*D SIN (AN) 

120 RP ( I ) =k (K) +DELV*DCGS ( AN) 

C 

C FCR THE UPPER FREE SURFACE 
DC L 30 1 = 55,62 

K= I NO ( I ,1) 

DELV= ( V ( I ) *"3- L . 0 ) * F AC U 
WRITE (6 , 121 ) K , D E L V 
XX ( I)=X(K)+DELV*DS1N(AN) 

130 RP ( I) =R (K)-DELV=CCOS ( AN ) 

122 FORMAT!/, » CORRECT IONS TO BE MACE AT EACH NODE ',/) 

121 FORMAT! 5X,‘ NOD E ' , I 5 , 5 X , « CORRE C T I CN ‘ , F 8 . 4 ) 
VC=V0-(V1-1.0)*V0 

c 

C RELOCATION OF THE NEW FREE SURFACE 
DC 140 1=9,62 
R = 3 

I F ( I .G E .55 ) M= 1 
K = I NO ( 1 , V ) 

X ( K ) = X X ( I ) 

R!K )=RR( I ) 

140 CONTINUE 

DC 160 I = 8C » 99 
I 1= INO (1,1) 

12 = 1 NO! 1,2) 

I 3=1 NO ( I ,3 ) 

IFd.GE.94) 12=13 

V3 = { 3 ( I 2 ) - B ! I 1 ) ) / ( X ( 1 2 ) -X ( I 1 ) ) 

W R I T t ! 6 , 1 6 1 ) 12, V3 

16 1 FORMAT ( 5X, »NODE= I5» 5X, ' VELOC= « ,F5.3) 

1 6 C CONTINUE 
C 

C EACH FIVE ITERATION CORRECT VO AND THE ANGLE BETA 
IF! INDEX. EC. 10) GO TO 101 
IF! INDEX. NE. 5) GO TO 100 
F ACU = 0 • 0 1 5 DO 
FACL = 0 . 00 IDO 
GO TO 100 
101 CONTINUE 

FACU=0. C01D0 
F ACL=0 . C 15C0 
C 

C CORRECTION CF THE ANGLE BETA 
DEE=O.03C0 
BET A = B ETA+DBE 
DC 19 1=55,59 
X C = 6 1- I 
K = I NO ( 1,1) 

X (K)=X( N2 ) -DX*XD 

19 R(K)=R(N2)-UTAN(5ETA)*(X(K)-X(N2)+DISP) 

ZR = R( 64 7 ) -R ( N2 ) 

X(617 )=X( 647)*0.25D0 
X (627 )=X (647) *0.500 
X(637)=X(647)*0. 75D0 
R(617)=R(N2)+ZR*0.25D0 
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R ( 627 ) =R(N2)+ZR*0.50C0 
R(637)=R<N2J+ZR*0. 75D0 
GC TO L GOO 
9999 CONTINUE 
STOP 
ENC 
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SLBRCUT ING BANDE C < A , F , NEC , M AXB » N'V E C , NROW , NG AN , NBEC I 
IMPLICIT REAL*8( A-r.O-Z ) 

DIMENSION A (NROW, MEAN) » F (NROW , NB EC ) 

LCCP=NEC-1 
CC 100 1 = 1, LOOP 

I F ( DAB S ( A ( I,1J J-L.E-30 ) 1000, 1000,10 
10 M = I + 1 

N=MINO( I + MAXb-1 , NEC ) 

DO 100 J = M,N 
L= J+2-M 

C=A(I»L)/A(I,1) 

DO AO M M= 1 , NVEC 
AO F { J , MM ) = F ( J , KM J -*D* F ( I , MM 1 
MM = MI MO (NAX8-L+1 , NEG-J + 1 ) 

CC 100 K=1,MM 
N N = L+K- 1 

IOC A(J,K) = A{ J,K)-D*A( I,NN) 

DC 199 1 = 1, NVEC 

199 F (NEQ, I ) = F (MEG, I )/A (NEC, 1 ) 

DO 20 0 1=2 , ME Q 

J=MEC-I +1 

l< = M I NO ( MEC-J + 1 , MAXE) 

DC 200 NM= 1 , NVEC 
DC 190 L = 2 , K 
N = J+L— 1 

19C F< J,MM) =F( J,MM)-A( J,L)*F(M,MM) 

200 F( J,MM)=F{ J,MM)/A(J, 1) 

RETURN 

10CC WRITE! 6, 1010) I 

1010 FORMAT ( 29ril ZERO DIAGONAL ELEMENT ROW IA) 

END 
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SIB ROUT INE DRVF IT < M C , L , X , Y , M , N , U , V , ANGL ) 



(Mi,Bi),(M2,B2),(M2,E2),(M4,B4), 
( X2 i PO ) , ( Y2 , QO ) , ( T 2 , G 1 1 
M1,M2»M3,M4 

(W2,Q2) ,(W3,Q3),(A1,P2),(B1,P3), 
( A 2 t D Z } ,(SN,R,Z) 



IMPLICIT REAL A 8( A-h,0-Z) 

C I MENS I CN X(L),Y(L) ,U(N) ,V(N) 

8 , AN GL ( N ) 

EQUIVALENCE 
I 

REALMS 
EQUIVALENCE 

I 

C PRELIMINARY PROCESSING 
1C M C C = M D 

MDML=MOC-l 
LC = L 

LMl=LO-l 
M C = Mi 

MM l = MO-l 
N C = N 

IF(MDO.LE.C) 

IF(MDO.GE.3) 

IF(LMl.LE.O) 

IF(MMl.LE.C) 

I FtNO.NE. LMl*MO+l ) 

GC TO ( II , 16) , MDC 
11 1 = 2 

I F ( X ( 1 ) -X { 2 ) ) 

12 DC 13 1 = 3, LO 

IF(X( 1-1 ) -X ( I) ) 

13 CONTINUE 
GC TO 18 

14 CC 15 1=3, LO 

I F ( X< I- 1 > -X ( I ) ) 

15 CONTINUE 
GC TO 18 

16 DC 17 1 = 2, LO 

I F ( X ( I - 1 ) . N E . X ( I J ) 

1 F ( Y ( 1-1! . E Q . Y ( I ) ) 

17 CONTINUE 

18 K = NO + MO 
I=LC+ 1 

DC 19 J= 1 , LO 
K=K— MO 
1 = 1-1 
U (K I =X ( I ) 

19 V(K)=Y(I ) 

RM = MO 

RM=1 .O/PM 
C MAIN CC-L OOP 
2 C K 5 = MO+ 1 

CC 80 1=1, LO 

C ROUTINES TO PICK UP NECESSARY X AND Y VALUES ARC 



GO TO 90 
GO TO 90 
GC TO 91 
GO TO 92 
GO TO 93 



12 ,95 ,14 
13,95,96 



96,55,15 



C-G TO 17 
GC TO 97 



TO ESTIMATE T F EM IF NECESSARY 



2 C 



21 



4 C 



I F ( I.C-T.l) 

X2 = U ( 1) 

Y3 = V ( 1) 

X4=U ( MO+1 ) 

Y 4= V ( M. 0+ 1 J 

A3=X4-X3 

E3=Y4-Y3 

IF ( MCM l . EQ . C ) 

I F ( LO .NE .2) 

A 4= A 3 
84= B3 

GC TO (33,32 ) , 
A2=A3+A3-A4 
A 1= A2 + A2- A3 
82=63+83-84 
E !=B2 + b2-b3 
GC TO (51,56) , 
X 2= X3 
Y2=Y 3 
X3=X4 



GO TO 40 



M3=B3/A3 
GO TO 41 



M CO 



M CO 
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GO TO 42 



V 3=Y 4 
X4 = X5 
Y4 = Y5 
A 1= A2 
Bl = 82 
A 2 = A 2 
62 = 83 
A 3 = A4 
63=84 

I F ( I .GE.LM1 ) 

41 K5=K5+MO 
X5 = U( K5) 

Y5 = V ( K5 ) 

A4 = X 5-X4 
64 = Y 5-Y4 

I F ( MCM1 .EC.O) 
GO TO 43 

42 IF(MDMl.NE.O) 
B4=B3+83-82 

42 IF(I.EQ.i) 

C-G TO (50,55 ) , 



M4=B4/A4 

A4=A3+A3-A2 



M CO 



GO TO 3 J. 



C NUMERICAL DIFFERENTIATION 

5 G T2=T 3 

51 W2=0A8S( R4-M3 ) 

N3=DAB$( K2-M1 ) 

S1m=W2 + W3 

IF(SW .NE.O.O) GO TO 52 

U 2= 0 . 5 
W 3 = 0 • 5 
SW= l . C 

52 T3= ( W2*M2+h3*M3 ) /SW 
IF(I-l) 80,00,60 

55 CCS2=CCS3 
S I N 2 = S I N 3 

5 t W2=CABS( A3*B4-A4~E3 ) 

R 3= D A B S ( A 1 8 2- A 2 ^ 6 1 ) 
IFIW2+W3.NE.0.01 GO TO 57 
K2=DSCRT <A3*A3+B3*B3 J 
N3=DSCiRT l A2*A2 + 6 2*B2 ) 

57 C 0 S 3 = ii 2 v A 2 + W3 * A 3 

S I iN3 = W 2 * B 2 + W3 ^ 8 3 
R=C0S3*C0S3+SIN3*SIN3 
IF(R. EG.O.O) GO TO 58 

R = t)SLR7 ( R) 

CCS3=CCS3/R 
S 1 N3 = S I N3/ R 

58 IF(I-l) 80,60,65 

C DETERMINATION OF THE COEFFICIENTS 

6 C G 2 = ( 2 .0*(M2-'T2 ) + M 2 - T 3 ) / A2 

G3= (-M2-M2+T2+T3 )/(A2*A2) 

GC TO 70 

£5 R=DS0RT(A2*A2+62*e2) 

F l = R ' C 0 S 2 

F2=3 . C*A2-R*(C0S2+C0S2+C0S3) 

P3=A2-P1-P2 

C1=R*SIN2 

C2=3. C*B2-R*( SI N2+S IN2+SIN31 
G3=B2-Q 1-Q2 
GC TO 75 

C COMPUTATION OF THE POLYNOMIALS 

7 C CZ=A2*RM 

2 = 0.0 

DO 71 J = 1 , MM 1 
„ K = K + 1 
Z=Z+DZ 
U ( K ) =P C+Z 

71 ' V(K) = G0 + Z*vGl + Z*(G2 + Z4Q3) ) 

GO TO 79 
75 1=0.0 

CO 76 J = 1 , MM 1 
K = K + 1 
Z=Z+RM 
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U( K ) = PC+Z*( Pl+Z*( P2 + Z*P3 ) ) 

76 V(K) = QC+Z*(Q1 + Z’MG2 + Z*Q2) ) 

79 K = K+ 1 

8C CONTINUE 
DC 82 1=1, 12 
82 ANGL ( I ) =0 . CDO 
NM I = NO- l 
CC 8 1 K = 1 2 i NM1 

81 ANGL( K ) =DATAN2( ( U ( K+ 1 ) -U ( K- 1 ) i » ( V ( K + 1 ) - V ( K - 1 ) ) > 

RETURN 
C ERRCR EXIT 



* 9 C 


NR ITE 


( 6,2090) 




GC TO 


99 


91 


WRITE 


(6,2091) 




GC TO 


9 9 


92 


W RITE 


(6,2092) 




GC TO 


99 


92 


WRITE 


( 6, 2C93) 




GC TO 


99 


95 


WRITE 


( 6,2095) 




GC TO 


98 


96 


WRITE 


(6,2096) 




GO TO 


98 


97 


WRITE 


(6,2097 ) 


98 


WRITE 


( 6,2093) 


99 


WRITE 


( 6,2099) 




RETURN 



I , X ( I ) ♦ Y ( I ) 
MDC,LO,MO,NO 



; FCRMAT STATEMENTS 
209C FCRMAT ( 1X/31H 

2091 FCRMAT ( 1X/22F 

2092 FCRMAT { IX/22H 
FCRMAT ( LX/25H 
FCRMAT { 1X/27H 
FORMAT ( 1X/33H 
FCRMAT ( 1X/33H 



2093 

2095 

2096 

2097 

2098 



FCRMAT ( 7h 



l 

2099 FCRMAT ( 7h 

1 2HM 

2 36H 
END 



**+ MD OUT OF PROPER RANGE./) 
*** L = 1 OR LESS./) 

*** M = 1 OR LESS./) 

*** IMPROPER N VALUE./) 

*** IDENTICAL X VALUES./) 

*** X VALUES OUT OF SEGUENCE./) 

*** IDENTICAL X AND 

I = , 14 ,ICX,6HX( I ) = , E 12 

IGX » 6HY ( I ) =,E12 
MC = , 1 9 , 8 X , 3 H L = , 15 , 8X, 

=iI5r8Xr3HK =,15/ 

ERROR DETECTED IN ROUTINE CRVFIT) 



Y VALLES./) 
3, 

3) 
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CCRNER NODAL POINTS AND THEIR COORDINATES 



MODE 


1 


X = 


-2.600 


R = 


0.0 


NODE 


3 


X = 


-3 . 6CC 


k = 


0.300 


NODE 




x= 


-3 . 600 


R = 


C.600 


NODE 


7 


x= 


-3.600 


R = 


0 .800 


MODE 


9 


x= 


-3 .600 


R = 


l .000 


NODE 


19 


X = 


-3.200 


R = 


0.0 


NODE 


21 


X = 


-3 .300 


R = 


0.300 


NODE 


22 


X = 


-3 . 3C0 


R = 


0.600 


NODE 


25 


x= 


-3.200 


R = 


0.800 


NODE 


27 


x = 


-3 . 300 


R = 


l .000 


NODE 


3 7 


x= 


-2.6C0 


R = 


0.0 


MODE 


39 


x= 


-2 .800 


R = 


0 .300 


NODE 


A 1 


x = 


-2 . 600 


R- 


0.600 


NODE 


A3 


X = 


-2 .600 


R = 


C.600 


MODE 


45 


X = 


-2.000 


R = 


1.000 


NODE 


c r; 
_•< 


X = 


-2.200 


R = 


0.0 


NODE 


57 


x = 


-2 .200 


n _ 
f\ — 


0.300 


NODE 


59 


X = 


-2 . 20 C 


R= 


0.600 


NODE 


El 


x= 


-2 .200 


R = 


C .800 


NODE: 


63 


x= 


-2 .300 


R = 


1.000 


NODE 


72 


X - 


-1 . 9CC 


R = 


0.0 


NODE 


75 


X = 


- 1 .90u 


R = 


0.300 


NODE 


77 


X = 


-1 .900 


R= 


C.600 


MODE 


79 


x= 


-1 .900 


R = 


0.600 


NCCE 


81 


X = 


- 1 .900 


R= 


1 .000 


NODE 


9 1 


X = 


-1 . 500 


R = 


0 .0 


NODE 


93 


X = 


-1.500 


R = 


C.3C0 


NODE 


95 


X = 


-1 .500 


R = 


0 .600 


NODE 


97 


x= 


-1 .500 


R = 


0.800 


NODE 


99 


x= 


-1.500 


R = 


1.000 


NCCE 


1 10 


X = 


-1 .3 50 


k= 


1 .000 


NODE 


1 1 1 


x= 


- 1 . 2C0 


R = 


0.0 


NCCE 


1 13 


X = 


- 1 .200 


R= 


0.300 


MODE 


115 


X = 


-1 .200 


K- 


0.600 


NODE 


117 


x= 


-1.200 


R = 


G.300 


NCCE 


1 19 


X = 


-I .200 


R = 


0 .97 0 


NODE 


12 1 


X = 


- 1 . 2 C 0 


R = 


1.000 


NODE 


125 


X = 


- 1 . 150 


R = 


1 .000 


NCCE 


130 


X = 


-1 .000 


R = 


0 .0 


NODE 


1 AC 


X = 


-I .CCC 


R = 


0.300 


NODE 


1-42 


X = 


- 1 .000 


R = 


0 .600 


NODE 


144 


X = 


-1 .000 


R = 


0.800 


NODE 


147 


x = 


- 1 . 100 


R = 


C . 9 7 0 


NCCE 


150 


X = 


-1 .100 


R = 


1 .000 


NODE 


153 


X = 


- 1 . C 3 C 


R = 


1.000 


NCCE 


156 


x= 


-1.000 


R = 


0.970 


NCCE 


158 


X = 


-1 .000 


R = 


1 .000 


NODE 


162 


X = 


-C .950 


R = 


1 .COO 


NCCE 


17 3 


X = 


-C .600 


R = 


0.0 


NCCE 


175 


X = 


-C . 80 C 


R = 


0 .300 


NOCE 


177 


x= 


— C .600 


R = 


C.600 


NCCE 


179 


X = 


-0 .600 


R = 


0.800 


NODE 


182 


X = 


- C . 9 C C 


R = 


0.955 


NOCE 


164 


x= 


-0 .800 


R = 


0.940 


NCCE 


189 


X = 


-0 .900 


R = 


1.000 


NOCE 


1 9 C 


x= 


-C .850 


R = 


1.000 


NCCE 


195 


X = 


-0 .76 5 


R = 


1 .005 


MODE 


2 C C 


X = 


-0 . 800 


R = 


1 .000 


NODE 


2 C 7 


x-- 


- C . 6 0 0 


R = 


0 .0 


NCCE 


2C9 


X = 


-C .600 


R- 


0 .300 


NODE 


2 1 1 


X = 


-C .600 


R = 


0.600 


NCCE 


2 12 


X = 


-0 .6' J 0 


R = 


0-814 
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NCDE 


217 


X = 


-0 .725 


R= 


0 .952 


NODE 


2 15 


x= 


-C «65L 


R = 


C .964 


NGCE 


224 


x= 


-0 .729 


R = 


UOU 


NCCE 


22 6 


X = 


-0.694 


R = 


]. . 0 1 9 


NCDE 


228 


x= 


- C . 66C 


R = 


1.028 


NCCE 


230 


X = 


-0 .630 


R = 


1 .037 


NCDE 


222 


X = 


-C .601 


R = 


1 .046 


NCCE 


244 


x= 


-C .400 


R = 


0.0 


NCCE 


24 6 


X = 


-0 .400 


R = 


0 .300 


NGCE 


2 4 £ 


X = 


- C . 4CC 


R = 


0.600 


NCCE 


25C 


X = 


-0 .472 


R= 


0.835 


NCCE 


253 


X = 


-0.587 


R = 


0 .931 


NOCE 


2 5 5 


X = 


-C . 522 


R = 


C.599 


NCCE 


26 T 


X = 


-C . 572 


R = 


1 .057 


NCCE 


20 3 


X = 


-C .544 


R- 


1 .069 


NGCE 


2 0 5 


x= 


-C.521 


R = 


1.030 


NCCE 


267 


X = 


-0 .493 


R = 


1 .092 


NGCE 


275 


X = 


C .C 


R= 


0.0 


NCCE 


2 £ 1 


x= 


-C . 100 


R = 


0.300 


NOCE 


233 


x= 


-0 . 160 


R = 


0 .600 


NCDE 


285 


x= 


-0.318 


R = 


C . 860 


NCCE 


286 


X = 


-0 .468 


R = 


1 .020 


NOCE 


25 C 


X = 


-C .414 


R = 


1.042 


NCDE 


256 


X = 


-C . 476 


R- 


1 .105 


NCCE 


258 


X = 


-0 .455 


R = 


1.120 


NODE 


3CC ' 


X =• 


-0.444 


R = 


1.128 


NCDE 
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